ANALYTICAL PROOF OF SPACE-TIME CHAOS IN 
GINZBURG-LANDAU EQUATIONS 



D. TURAEV AND S. ZELIK 

Abstract. We prove that the attractor of the ID quintic complex Ginzburg-Landau 
equation with a broken phase symmetry has strictly positive space-time entropy for an 
open set of parameter values. The result is obtained by studying chaotic oscillations in 
grids of weakly interacting solitons in a class of Ginzburg-Landau type equations. We 
provide an analytic proof for the existence of two-soliton configurations with chaotic 
temporal behavior, and construct solutions which are closed to a grid of such chaotic 
soliton pairs, with every pair in the grid well spatially separated from the neighboring 
ones for all time. The temporal evolution of the well-separated multi-soliton structures 
is described by a weakly coupled lattice dynamical system (LDS) for the coordinates 
and phases of the solitons. We develop a version of normal hyperbolicity theory for 
the weakly coupled LDS's with continuous time and establish for them the existence of 
space-time chaotic patterns similar to the Sinai-Bunimovich chaos in discrete-time LDS's. 
While the LDS part of the theory may be of independent interest, the main difficulty 
addressed in the paper concerns with lifting the space-time chaotic solutions of the LDS 
back to the initial PDE. The equations we consider here are space-time autonomous, 
i.e. we impose no spatial or temporal modulation which could prevent the individual 
solitons in the grid from drifting towards each other and destroying the well-separated 
grid structure in a finite time. We however manage to show that the set of space-time 
chaotic solutions for which the random soliton drift is arrested is large enough, so the 
corresponding space-time entropy is strictly positive. 



1. Introduction 

We demonstrate that if an evolutionary system of partial differential equations (PDE) 
in unbounded domain has a solution localized in space and chaotic in time, then one 
should expect both temporal and spatial chaotic behavior in the system. Namely, one 
may observe a formation of non-trivial spatial patterns that evolve in an irregular fashion 
with time, and the corresponding space-time entropy [CoEc99j IZel07] is strictly positive. 
In other words, the number of solutions which are essentially different from each other on 
a finite space-time window grows exponentially with the window volume. 

As a tool for finding the spatially-localized, temporally-chaotic solutions one may try, as 
we do it here, to look for special types of both spatially and temporally localized solutions. 
Thus, like Shilnikov homoclinic loop and Lorenz butterfly serve as a criterion for chaos 
formation in systems of ODE's |Sh i65l IS~hi70l IShSTl ITuSh] . the existence of the Shilnikov 
homoclinic loop in the dynamical system generated by the PDE on the space of spatially 
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localized solutions implies the space-time chaos in the extended system that corresponds 
to uniformly bounded solutions of the same PDE. 

We do not prove this principle in full generality here. Instead, we decided to show 
how it works for a class of Ginzburg-Landau equations with a broken phase symmetry. 
The main motivation for such approach is that despite a huge amount of numerical and 
experimental data on different types of space-time irregular behavior in various systems, 
there are very few rigorous mathematical results on this topic and mathematically relevant 
models describing these phenomena. Therefore, we made an effort of providing a free from 
numerics, completely analytic proof of the existence of space-time chaos in an important 
equation of mathematical physics. 

The basic mathematical model for the space-time chaotic behavior is the so-called Sinai- 
Bunimovich chaos in discrete lattice dynamics, see [AfrFeOOt IBuSi88t IPeSi88l IPeSi91j . 
This model consists of a Z n -grid of discrete-time chaotic oscillators coupled by a weak 
interaction. The single chaotic oscillator of this grid is described, say, by the Bernoulli 
scheme Ai 1 '■— {0, 1} Z , so the uncoupled system naturally has an infinite-dimensional hy- 
perbolic set homeomorphic to the multi-dimensional Bernoulli scheme A4 n+1 := {0, 
(Ai l ) zn . The temporal evolution operator is then conjugate to the shift in Ai n+1 along 
the first coordinate and the other n coordinate shifts are associated with the spatial trans- 
lations on the grid. Due to the structural stability of hyperbolic sets, the above structure 
survives under a sufficiently weak coupling. Thus, in this model, the space-time chaos is 
described by the multi-dimensional Bernoulli scheme M. n+1 . 

Importantly, the space-time entropy in the Sinai-Bunimovich model is strictly posi- 
tive. We know from the general theory of dissipative systems in unbounded domains 
(see e.g. |CoEc99t IMirZe08l IZel03t IZel04t IZel07j ) that under some reasonable dissipativ- 
ity assumptions this entropy is finite for systems of evolutionary PDE's, therefore the 
Sinai-Bunimovich model carries "enough complexity" to be able to capture certain basic 
features of spatio-temporal chaos in systems of various nature. In particular, it is well 
established by now (see e.g. |Pom86t IMan09t IBaT05j ) that the transition from regular 
to chaotic space-time behavior often happens via the emergence of well spatially sepa- 
rated and long living "turbulent spots" . As the interaction between such spots seems to 
be weak, the Sinai-Bunimovich chaos paradigm can be relevant for the analysis of these 
near-threshold phenomena. 

Yet, a direct application of the Sinai-Bunimovich construction to systems with contin- 
uous time and space is not possible, in general. Even the existence of one PDE which 
possesses an infinite-dimensional Bernoulli scheme was a long-standing open problem. 
The first examples of such PDEs (in the class of react ion- diffusion systems), have been re- 
cently constructed in |MieZe07] . The method used in that paper is based on a strong and 
explicit spatio-temporal modulation of the equation right-hand sides, which effectively 
transforms the systems into a discrete-time lattice dynamical system. The disadvantage 
is that very special (and artificial) nonlinear interaction functions emerge in the result, 
which are far from the usual nonlinearities arising in physics models. 

A different approach to the problem is suggested in |MieZe09"] , where a theory of weak 
interaction of dissipative solitons was developed and, as an application, a space-time 
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chaotic pattern has been constructed for the perturbed ID Swift-Hohenberg equation 



(1.1) d t u + (d 2 x + l) 2 u + (3 2 u + f(u) = fih(t,x,fj), f(u)=u 3 + nu 2 , n < 1. 



Here /i/z(t, x, fi) is a space-time periodic forcing. Its exact form is quite non-trivial, how- 
ever the amplitude /i can be taken arbitrarily small. The idea is to create a spatially 
localized spot of chaotic temporal behavior, and to build then a grid of such spots, well 
separated in space. The spots are pinned down to the prescribed locations at the grid 
points by spatial oscillations in the forcing [ih. If the spots stay sufficiently far apart, 
their interaction is small, so a small amplitude forcing occurs to be sufficient to sustain the 
grid for all times (the wave length of the forcing has, however, to grow as the amplitude 
decreases) . 

Equation (jl.ip at fi = 0, like many other important equations, does have a spatially lo- 
calized solution, a soliton, u = U(x) with exponentially decaying tails. One may therefore 
look, at all small fi, for multi- soliton solutions in the form 



where £j(t) is the position of the j-th soliton; the well-separation condition reads as 
L := ir&j-Lk — £fc|| ^ 1- Due to the "tail" interaction and the small forcing, the solitons' 
positions £j(t) may move slowly, and this motion is described by a lattice dynamical 
system (LDS), see |MieZe09] for details. The obtained LDS is not in the form one needs 
for establishing the Sinai-Bunimovich chaos (a grid of chaotic maps with weak coupling), 
since the individual solitons u = U(x) are equilibria at fi = and do not have their own 
(chaotic) dynamics. However, as it is shown in |MieZe09] . a pair of weakly interacting 
solitons in the ID Swift-Hohenberg equation can be forced to oscillate chaotically in 
time by an appropriate choice of the time-periodic perturbation jj,h{t,x,n). For a well- 
separated grid of such soliton pairs, one obtains a time-periodic LDS, and the period 
map for this system is the sought discrete lattice of weakly coupled chaotic maps, i.e the 
space-time chaos is established. 

The scope of |MieZe09] is much more general than the Swift-Hohenberg equation: by 
developing the center manifold approach proposed in [San02j, the paper derives the LDS 
that governs the evolution of weakly coupled multi-soliton configurations for a large class 
of systems of evolutionary PDE's. It also proposes a method for constructing spatially 
localized and temporally chaotic solutions which are obtained as a system of finitely many 
weakly coupled stationary solitons. Note that, although spatially localized solutions with 
non-trivial temporal dynamics have been observed numerically and experimentally in 
various physical systems (see e.g. |AkSCT0"T| IBiWe02[ ITuVZ07j and references therein), 
the direct analytic detection and study of such solutions is obviously a very difficult task. 
However, when a finite system of well-separated solitons is considered, the descriprion 
provided by [MieZc09j for the evolution of such object is often a low- dimensional system 
of ODE's which can exhibit a chaotic dynamics |TuVZ07] and can be studied analytically, 
so the chaotic temporal behavior of such localized patterns can be rigorously proven. 

In the present paper we show how a space-time chaotic lattice can be built out of these 
chaotic multi-soliton systems in the case where no spatial nor temporal modulation is 
imposed. Two problems immediately appear in this setting: 
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1. with no external forcing, the LDS which describes the multi-soliton dynamics is an 
autonomous system with continuous time, so the Sinai-Bunimovich chaos construction 
(for which the discreteness of time is very essential) is not applicable; 

2. with no spatial modulation, there is no pinning mechanism which would keep the 
solitons eternally close to any given spatial grid, therefore the infinite-time validity of the 
LDS description is no longer guaranteed. 

We resolve here both the issues. As an application, we consider the ID quintic complex 
Ginzburg-Landau equation with slightly broken phase symmetry: 

(1.2) d t u = (1 + ip)dlu - (1 + i5)u + (z + p)\u\ 2 u - (e 1 + ie 2 )\u\ A u + P, 

where /3, 5, p, Eip, p are some real parameters, and /i < 1. We mention that, in contrast to 
the previous on the Swift-Hohenberg equation, we do not have here any artificial functions, 
and the only freedom we have is the choice of the numeric parameters. Note also that 
the Ginzburg-Landau equation serves as a normal form near an onset of instability, i.e. it 
very often appears in applications as a modulation equation for various more complicated 
problems. The phase symmetry in the modulation equation appears as an artefact of 
closeness to the instability threshold, so if there is no such symmetry in the original 
problem, then the effects of small symmetry breaking also need to be considered, see 
[Mie02j and references therein. While we introduce only the simplest symmetry breaking 
term ("+//") in (11.21) . the general case is also covered by the theory (see Section [2]). 
The main result of the paper is the following theorem (Section [3]). 

Theorem 1.1. There exists an open set of parameters (/?, S, p, Ei, 62, p) such that equation 
(II. 2p possesses a global attractor A (say, in the phase space Lf(M.)) with strictly positive 
space-time entropy 

h,- t (A) > 0. 

Equation (ll.2p at p. = has the additional phase symmetry u — > e l ^u. Therefore, for 
each stationary soliton u = V(x) of this equation, u = e^V(x) is also a stationary soliton. 
Therefore, the multi-soliton configurations are given by 

u(t,x) = ^2e iMt) V(x - + "small corrections", 

i 

where £j and <pj are the coordinate and phase of the j-th soliton. For a soliton pair with 
the states (£1, 0i) an d (£2, 02), the evolution is governed, to the leading order with respect 
to the distance \£ 2 — by the following system of ODE's: 
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see [VlKROll ITuVZ07l IMieZe09j . Here r is a scaled slow time, R := (£ 2 - £i)/2, $ : = 
(pi — 02, ^ : = (0i + 02)/2 and a, b, u>, 6^2, c, v and Q are parameters whose exact values 
depend on the values of the original parameters of (II. 2p (see the corresponding expressions, 
as well as asymptotic expansions near the exactly solvable nonlinear Schrodinger equation, 
in Sections I2H3I) . While the variables R, $ and \l/ can be treated as the "internal variables" 
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of the two-soliton pattern, the variable p := (£1 + ^)/2 marks the spatial position of the 
soliton pair. To the leading order, it is governed by the equation 

(1.4) -f-p = -e~ aR cosicoR + 6 X ) sin($) := g(R, $). 

dr 2 

A numerical study of system (II. 31) undertaken in |TuVZ07] revealed various chaotic 
regimes for different parameter values. In order to provide an analytic proof (see Lemma 
12. 3p of the chaotic behavior (i.e. the existence of a nontrivial hyperbolic invariant set) in 
this system for an open set of parameter values, we have found a degenerate equilibrium 
of that system with 3 zero eigenvalues. The normal form equation for this bifurcation 
(see |ArCST85] ) leads to the following 3rd order equation: 

Y'» = 1 - Y 2 + EY', 

where E (a certain combination of the parameters of the original system) can take any 
real values. In [IbRo05], the existence of a Shilnikov homoclinic loop for this system was 
proven at certain E values, which implies [Shi65t IShi70] chaos for some interval of the E 
values and, hence, for an open set of parameter values for system (II. 3p . 

Chaotic solutions of system (II. 3p correspond to a chaotically oscillating soliton pair, 
which is a temporally chaotic and spatially localised solution, by construction. After that, 
according to the program described above, we build a well spatially separated lattice of 
such time-chaotic solitons. The center manifold reduction theorem proved in |MieZe09] 
ensures that the evolution of this lattice is governed by a system of infinitely many weakly 
coupled copies of the ODE's fO]) . fOj) . 

Even when every individual ODE-subsystem in the continuous time LDS is hyperbolic, 
the LDS itself is not hyperbolic (this is a principal difference with the Sinai-Bunimovich 
chaos in the discrete-time LDS's where the countable product of hyperbolic sets for the 
individual maps is hyperbolic again). Each constituent ODE contributes a neutral direc- 
tion corresponding to the time shift, so for the linearized flow of the continuous time LDS 
we have infinitely many neutral directions. Therefore, after a weak coupling is switched 
on, the dynamics is not preserved (the LDS can hardly be topologically conjugate to the 
uncoupled one). Still, the invariant manifold theorem of Section H] shows that if, given any 
orbit of the uncoupled LDS, we consider the family of all orbits obtained by all possible 
time-reparametrizations in each of the constituent ODE's, then this family continues in a 
unique way as an invariant manifold of the weakly coupled LDS. This fact allows to show 
the strict positivity of space-time topological entropy for the countable systems of weakly 
coupled chaotic oscillators with continuous time. 

In fact, results of Section H] cover LDS's of a more general type. The problem we 
have to deal with is that, although system (ll.3p for the internal variables (R, $, \l/) of the 
chaotic soliton does have a uniformly hyperbolic set, the full system describing the motion 
of the chaotic soliton includes equation (jl.4p for the soliton position p, and is clearly 
non-hyperbolic (so we have to consider the LDS's built of partially-hyperbolic individual 
ODE's). The neutral directions appear because the right-hand sides of (ll.3p .( ]l~4"]) are 
p-independent, which is a mere consequence of the translational symmetry of the PDE 
under consideration, i.e. their presence is an inherent property of the soliton-interaction 
equations in systems without a spatial modulation. 
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Since the internal variables of the soliton change chaotically with time, the soliton 
position p(t) performs, essentially, an unbounded random walk (as an integral of a chaotic 
input, see (11.4p ). When the chaotic solitons are well spatially separated, the contribution 
of the neighboring solitons to the p-equation is small, so we have essentially independent 
random walks for each of the chaotic solitons in the lattice. This makes it impossible 
for us to ensure that the well-separation condition is fulfilled for all times and all initial 
multi-soliton configurations. We, in fact, believe that the majority of these configurations 
do break up this condition in a finite time, so the corresponding solutions cannot be 
completely described by the weak soliton interaction paradigm. 

However, the weak soliton interaction theory of |MieZe09] is the only tool we have here 
for the analysis of dynamics of the multi-soliton patterns. As we are unable to control the 
soliton's random walk, we devise a method of keeping track of those configurations for 
which the well-separation condition holds eternally (i.e. the LDS description is applicable). 
This method allows us to verify (Section [6]) that the number of such solutions is large 
enough to ensure the positivity of the space-time entropy. It is worth to emphasize that, 
instead of fighting with the random walks, our method exploits them in a crucial way. 

Roughly speaking, assume that the hyperbolic set for (II. 3p contains two periodic orbits 

l\ and T 2 and a number of heteroclinics which connect them. Assume that, according 

to equation (11 .4p . the soliton pair moves to the left if (R, belongs to I\ and to the 

right when it belong to T 2 . The direction of this motion is determined by the sign of 

1 f Tj ... 
bj ■= — g(R J (t), &(t)) dt, where 7} is the period of Tj = (R J , $ J , ^). So, we require 
T j Jo 

bib 2 < (in fact, only b\ ^ b 2 is enough, as we show). Then, our orbit selection method 
works as follows: assume that initially the j-th soliton is in the interval [L~ , Z/t] with 
Li- — Lj large enough; then until it remains in that interval, we allow the internal state 
(R, $, \I/) of the soliton to jump randomly between l\ and T 2 along the heteroclinic orbits 
(thus we gain the complexity which is enough to have the positive entropy); however, 
when the soliton reaches the bound (say, Lj), we stop allowing jumps and consider only 
orbits that stay near T 1 until the soliton position Pj(t) arrives close to (Lf + Lj)/2 (when 
the bound LJ is achieved, the orbit must stay near T 2 ); after Pj(t) is driven to the middle 
of the interval, the random motion is allowed again, and so on. 

We proved in Section that the above described procedure can be implemented simul- 
taneously for all chaotic solitons on the grid, and it allows indeed for a selection of a set 
of spatially non-walking solitons with positive space-time entropy. In order to do this, 
we need a further development of the theorem on normally-hyperbolic manifolds in the 
countable product of partially hyperbolic sets which is proved in Section HJ namely we 
prove certain "asymptotic phase" results in Section [5j 

As the above discussion shows, the theory we build is readily applicable to any dissi- 
pative PDE for which the weak soliton interaction system for some finite multi-soliton 
configuration exhibits a chaotic dynamics. In analogy to the finite-dimensional case, we 
are now able to analyze localized structures and effectively use them for the understanding 
of space-time dynamics generated by PDEs. 

We are grateful to A.Mielke and A.Vladimirov for useful discussions, and to WIAS 
(Berlin) and BenGurion University for the hospitality. 
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2. Space- time chaos in complex Ginzburg-Landau equation with broken 

phase symmetry 

Consider the one-dimensional complex Ginzburg-Landau equation 

(2.1) d t u = (1 + i(3)d 2 x ii - (1 + i8)u - uH(\u\ 2 ) + (iG{u) 

where u = U\ + iu^ is an unknown function of x G R and t G R, the function H : R — > C 
is smooth, i?(0) = 0, and parameters /3, 5 are real; the symmetry-breaking parameter fi 
is assumed to be small, and the function G is smooth. 

Let K, be a set of solutions of ( 12. ip which are defined and uniformly bounded for all 
(t, x) G R 2 (under certain standard dissipativity assumptions, equation (12. ip will have 
global attractor; in this case one can choose as the set K. the set of all solutions that lie in 
the attractor, see more after Theorem 12. ip . The complexity of spatio-temporal behavior 
of the solutions can be characterized by the space-time topological entropy defined as 

(2.2) h s AK,) = lim limsup — hJ)C\ ul ^„, 

K J K J £-►() (Tifl) _>£ 4772 V l|t|<T,M<^ 

where /CL, stands for the set of functions from K, restricted on the space-time 

i\t\<l ,\x\<H 1 

window {\t\ < T, \x\ < R}, and h e denotes the Kolmogorov e-entropy of this set, i.e. the 
logarithm of a minimal number of e-balls in the space L°°([—T,T] x [— R,R]) which are 
necessary to cover Q the set; see |KoTi61j . It is well-known (see, e.g., |CoEc99l IMirZe08l 
IZel04(. [Zel07J ) that the space-time topological entropy h s _ t (JC) is well-defined and finite in 
our case. Thus, if h s _ t ()C) is strictly positive for some set /C, then the number of various 
spatio-temporal patterns that are supported by the equation grows exponentially with the 
volume of the space-time window. 

In our construction of spatio-temporal chaos we assume that the nonlinearity H is such 
that for some (3 = /3 and S = So the Ginzburg-Landau equation ( 12.1 ft possesses at \i = 
a stationary, spatially localized solution u = U(x): 

(2.3) (1 + ip )d 2 x U - (1 + iS )U - UH{\U\ 2 ) = 0; 

for the existence results see [AfMi99t [AfMiOl] and references therein, and Theorem 13.11 

Equation (12. ip is invariant with respect to spatial translations x — > x — £ and, at fi = 0, 
with respect to phase shifts u — > e^u. So, along with the given soliton U(x), equation 
(12. ip possesses at fi — a family of stationary solitons: 

(2.4) u = U^(x) := e^U(x - £), (6 0) G R 1 x S 1 . 

Because of the symmetry with respect to x — > —x, along with the soliton u = U(x), 
equation ( 12. 3 p also has a localized solution u = U(—x). Equation (12. 3p is an ODE with 
4-dimensional phase space. A localized solution corresponds to a homoclinic intersection 
of the stable and unstable manifolds of the zero equilibrium of this system. Since these 
manifolds are 2-dimensional and family ( 12 .4p is 2-parametric, all the localized solutions of 
( 12. 3 p are contained in family ( 12 .4p . Thus, U(—x) = e l<t>0 U(x — £ ) for some 0o>£o> which 



it follows in a standard way from the parabolic regularity, that K.\ i t , <T \ x \<r i s compact 
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immediately implies that U^/2,o{~ x) — ^U^/ 2} o(x). In other words, we may from the very 
beginning assume that our soliton is chosen such that it is either symmetric: 

(2.5) U(-x) = U(x), 

or antisymmetric (U(—x) = —U(x)). In this paper we consider the symmetric case, i.e. 
we assume that (12.51) holds (in the antisymmetric case the soliton interaction equations 
are different; however one can show that a small perturbation of an equation with anti- 
symmetric soliton creates symmetric solitons - cf. [AfMi99, AfMiOl], so the results of our 
paper can be applied in this way). 

Since every function in f)2.4p is a stationary solution of (12.31) at fi = 0, it follows that 



the functions (pi := —d^U^o ^ =0 = 9 X U and ip 2 := d^Uo^^ = iU belong to the kernel of 
the linearization C v of (12. 3p at U: £^1,2 = 0, where 

(2.6) C v <p := (1 + ifo)%<p - (1 + i5o)<p - H(\U\ 2 )^ ' \U\ 2 H*{\U\ 2 )<p - U 2 H'{\U\ 2 )^ 

(<p is a complex conjugate to (p). Thus, zero is a double eigenvalue of £jj- 

We assume that the soliton U is non-degenerate in the sense that the rest of the spectrum 
of £u is bounded away from the imaginary axis; e.g. the algebraic multiplicity of the zero 
eigenvalue is two (note that since U(x) — > as x — > ±oo, the operator L v is a compact 
perturbation of the operator <p \-t (1 + i(3 )d 2 ip — (1 + iS )(p, so the essential spectrum is 
bounded away from the imaginary axis; however, one should check that the eigenvalues 
stay away from the imaginary axis as well). 

Under the non-degeneracy assumption, the conjugate operator C\j, which we define as 

(2.7) CU := (1 + iPo)S^ - (1 + *W - H(W?)^ - \U\ 2 H'(\U\ 2 )4> - U 2 H\\U\ 2 )4>, 

also has a two-dimensional kernel. The corresponding pair of adjoint eigenfunctions ipi 
and ip2 can be chosen such that 

/+oo 
(Pi(x)ipj(x)dx = Sij, ipi(-x) = -ipi(x), ip 2 (-x) = ip 2 (x). 
-oo 

As x — > ±oo, the functions U, (pi, ipi decay exponentially, with the rate A given by 

(2.9) ReA = -a<0, ImA = w, (-a + iu) 2 (l + i/3 ) = (1 + iS ), 
see [AtMi99L IMieZe09j for details. Thus, we have 

(2.10) ~re ( - a+iu)|x| , ~ se^ a+iui)lxl sign(x), ip 2 ~ qe ( - a+iuj)lxl as |z| -> oo, 
where r, s, q are some non-zero complex constants. We introduce the notation 

(2.11) ae idl := 4isr(l + i/3 )X, be id2 := Atqr(l + i/3 )X, 9 := 6 2 - B x . 
Denote 

/+oo 
e-^il) 2 {x)G{e i4> U{x))dx, 
-oo 

where G(u) is the symmetry-breaking term in ( 12.11) . Since F(<f>) is periodic, the equation 
(2.13) F'(<f>* + j) + F'(f - j) = 
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always has solutions. We assume that there is a solution <p* such that 

(2.14) c := 2F'(0* + ^ 0, 

(2.15) F "(0* + ^) + F "( *_!^o. 
Conditions fl2.13p - fl2.15p define the constant 0*. Denote also 

(2.16) 7 := V(0* + ^)-F(0*-^)]. 

c 4 4 

In the basic case G(u) = 1, we have F(<f) = ccos(</> — Q, where ce^ = if)2(x)dx. It 
is easy to see that 0* = £, c = — ca/2 and 7 = in this case, and that both conditions 
f l2~14j) and fl2TT5|) are fulfilled provided ^^2{x)dx ^ 0. 

Theorem 2.1. Lei, along with (I2.14p . f)2.15p . the following conditions be satisfied for a 
non- degenerate, symmetric stationary soliton U(x): 

(2.17) a ^ 0, b j£ 0, u ^ 0, u ^ 2 7 ct, cos 9^0, a sin 9 + w cos 9 ^ 0, 



a 



2 



(2.18) 4a;-(cos6' + 27sin6») < |l + 2 7 -(a cos 6» + u sin 6» 

Then, arbitrarily close to fi = and <5 = <5o i/iere exist an interval of values of fi and an 
interval of values of 5 such that the corresponding equation (12. ip has a uniformly bounded 
set of globally defined solutions with strictly positive space-time entropy. 

Proof. Each of the solutions of equation (12. ip that belong to the large (of positive entropy) 
set we are going to construct can be viewed as a slowly evolving multi-soliton configu- 
ration. Namely, we choose a sufficiently large L and consider solutions u(x,t) which for 
every t G R stay close, in the space C&(R) of bounded continuous functions of x, to the 
multi-soliton manifold defined as the set of all functions u(x) of the form 

(2.19) u(x) = u m :=J2 U ^ : " 

jez 

where m := is any sequence such that 

(2.20) inf(£ i+ i-&) >2L. 



For sufficiently large L, the multi-soliton manifold is indeed an infinite-dimensional sub- 
manifold of C&(R) which is parameterized by the sequences m := {^j, <j)j} of the soliton po- 
sitions and phases (see [MieZe09] ). The boundary dM.L is given by infj £ z(Cj+i = 2L. 

We will seek for solutions of equation (12. ip in the form u(t) := w m ( t ) + w(t) where m(t) 
is a slow trajectory in and w(t) is a small corrector. Recall a result from [MieZe09j. 



Theorem 2.2. For all L large enough there exists a C k -map S : Ml — > Cb(R) such that 

\C k (} 



(2-21) IISIIcfcfM^GVR)) < Ce aL 
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(where a > is defined by ( 12.91) ) and that the manifold S := {u = u m + S(w m ), m G M L } 
is invariant with respect to equation (12. ip . Namely, there exists a C k -vector field T on 
Ml such that given any solution of 

(2.22) |m(t) = ^(m(t)) 

defined on a time interval t G (t_,i+) ; the function 

(2.23) u m (t) +S(«m(t)), te(L,t + ), 

solves equation (12.11) . 

Moreover, system (12.221) /ias £/ie following form: 

— & = 2Re[ s r (1 + iB a )\ { e *&+i-&)+*(fc+i-*i) _ e Afo-fc-i)+i(*j-i-fc)}i + . 

( 2 - 24 ) 1^. = _ 2 Re[ g r (1 + i(5 )\ [ e x ^~^ +i ^-^ + e^-fc-O-Hte-i-fc) }] 

+/iF(0 J ) - (5 - (Jo) + • • • , 

where a and co are the same as in (12. 9p . the constants r,s,q are defined by (12.101) . the 
function F is defined by (12 . 1 2[) . and the dots stand for terms which are 0(e~ 3aL + /i 2 + 
(5 — So) 2 ) in C k (M.L, M.) -metric, uniformly for all j G Z. 

The complete proof of formulas (12.241) occupies a substantial part of |MieZe09"] . For 
reader convenience, we provide a brief heuristic derivation of the equations in Appendix. 

According to Theorem 12.21 the evolution of well-separated multi-soliton configurations 
in the driven Ginzburg-Landau equation is governed by system (12.241) . Therefore, in order 
to prove the positivity of space-time entropy in equation ( 12.11) . it is enough to find a large 
set of solutions of system (I2.24H which satisfy the separation condition ( I2.20p l^ 

The corresponding theory for a class of lattice dynamical systems which includes system 
(I2.24p is built in Sections HH6J In particular, Theorem 16.11 gives a general result on the 
existence of a set K. of non-walking trajectories of a lattice dynamical system such that 
h s _ t {JC) > 0. In what follows we will show that a certain subsystem of (I2.24p indeed 
satisfies conditions of Theorem 16.11 

For any v and any sufficiently large L we may define a sequence L n , n G Z, as follows: 

(2.25) L = ve~ aL t, L 2n+1 = L 2n + 4L, L 2n+2 = L 2n+1 + 2L. 
We will look for pulse configurations which satisfy £j(t) = Lj + r)j(t) where 

(2.26) \r]j{t)\ <C, t G M, j G Z, 

where the constant C is independent of L, j and t. In other words, we have a grid of 
weakly interacting pulse pairs with the distance between the pulses in the pair of order 
2L and the distance between pairs of order 4L. Assumption (12.261) then means that we 
should ensure that this structure is preserved for all t although a uniform spatial drift of 
the whole grid is allowed (ve~ aL is the velocity of the drift). 



2 recall that system (|2.24j) is defined on the manifold whose boundary is given by (|2.20j) : outside 
this boundary the reduction to the invariant manifold S may fail - the so-called strong soliton interaction, 
soliton collisions, etc., may take place 
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Further, we introduce the scaling r := te~ , O := (S — 5 )e , v := fie . We 
will consider a region of bounded fl and i/, which corresponds to fi and 5 — 5q of order 
0(e~ aL ). We also assume L = — , n G N. Equations (I2.24[) recast as follows (see f)2.1ip ): 



(2.27) 



d 






= v - 


d 


v + 


dr^' +2 = 


d * 


b 


dr"^ 1 


~ 2 


d * 


6 




~ 2 



e -««i sin(a;i2j - ^ + X ) + 0( e - Qi ), 



2 

—aRn 



-aL\ 



sin(uRj - + 2 ) + vF{(j) 2j+ i) - Q + 0(e 
sm(cuRj + $j + 2 ) + vF(<f) 2j+2 ) - tt + 0(e 



-aL\ 



-aL\ 



where we denote Rj := 772.7+2 — V2j+i, $j '■= — (4>2j+2 — <f>2j+i)- As we see, only interaction 
inside the soliton pairs gives a contribution into the leading terms of equations ( 12.271) : 
since the distance between pairs is, in our configuration, of order 4L, the leading term 
for the interaction between solitons from different pairs will be of order 0(e~ 4aL ) in the 
non-rescaled time t, so after the time rescaling it is of order 0(e~ 2aL ), i.e. it is absorbed 
in the 0(e" aL )-terms in (|2T27|) . 

Let us rewrite the system in the coordinates Rj, tyj := (4>2j+i + 02j+2)/2, and 
Pj '■= (V2j+i +V2j+2)/2 (i-e. Pj is the center of the soliton pair, Rj is the distance between 
the solitons in the pair, and ^ describe the soliton phases). We obtain 



(2.28) 
(2.29) 



-^-Pj = v + aR] cos (a; Rj + 6>i) sin($ J ) 



0(e~ aL ) 



dRi 



ae- aR i sm(uRj + 6 1 )cos($ j ) + 0(e' aL ), 



dr 
d<£> ■ 



be~ aR >cos(cuRj + 9 2 )sm($j) + v 



dr 



2 



v 



k dr 



sin (a;i^ + 2 ) cos ( $■,■)+ - 



n*j+-£)-F{*j -2 



-fi+0(e 



At large L this system is a lattice dynamical system of form f)4.9p : at L — +oo the 
subsystems that correspond to different j are independent and identical, and the equations 
for variables fjj := (Rj,$j,tyj) (equations (I2.29P ) are independent of the p-equation 
(E2H1) 



Therefore, in order to prove Theorem 12.11 it is enough to check that system 
(I2.28p .f l2.29p satisfies conditions of Theorem 16.11 at some v. According to that theorem, 
we will then obtain, for all sufficiently large L, the existence of a set K, of solutions of 
system fl2.28p . fl2.29p which has a positive space-time entropy and is uniformly bounded by 
a constant independent of L (i.e. condition f !2.26|) is fulfilled - this, in turn, ensures that 
the separation condition fl2.20p holds, with somewhat smaller L, for all the solutions from 
K,). Now, lifting the set K by formula f)2.23p . we obtain a uniformly bounded set /C of 
globally defined solutions of the perturbed Ginzburg-Landau equation, and the positivity 
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of the space-time entropy of the set /C follows from the smallness of § and the positivity 
of the space-time entropy of K. 

Thus, to finish the proof we need the following 

Lemma 2.3. Assume (I2.13p - (12.18p . Then there exists an open region of values of v and 
Q for which the system 



d 

(2.30) 



ae 



-aR 



sm(uR + 9i) cos($), 



$ = be- aR cos(uR + 9 2 ) sin($) + v 



dr 
d 

dr 



e- aR sm(uR + 9 2 )cos(<$>) + 



v 



F(¥+-)+F(¥-- 



■n 



behaves chaotically, i.e. it has a basic (= non-trivial, uniformly-hyperbolic, compact, 
locally-maximal, transitive, invariant) set A. Moreover, in A one can find two periodic 
orbits, y = y-(r) and y = y + (r), of periods T_ and T + , respectively, such that 



(2.31) 



1 

fZ 

aR , 



g{y-{i~))dr ^ 



where g(y) :— v + |e cos(uR + 9\) sin($) . 

One may check that condition (I2.3ip implies that 



(2.32) 



g{y-{i~))dr 



g(y + (T))dr < 



for an appropriately chosen v. Hence, the lemma indeed establishes the required fulfilment 
of conditions of Theorem 16.11 chaotic system f l2.30p coincides with the y-subsystem ( I2.29P 
at L — +oo (for every j), and condition f !2.32p coincides with condition (I6.22p (the function 
g is the right-hand side of the p-equation f)2.28p ). So, it remains to prove the lemma. 

We note that numerically the existence of chaos in system (I2.30p with F = ccos^> as 
well as different scenarios of its emergence for various parameter values were established in 
|TuVZ07] . In our analytic proof of chaotic behavior we use one of the scenarios mentioned 
in |TuVZ07] . Namely, we find an equilibrium of system (I2.30p with 3 zero characteristic 
eigenvalues. It is known |ArCST85| IIbRo05] that bifurcations of such equilibrium lead to 
a Shilnikov saddle-focus homoclinic loop, hence to chaos. 

Proof of Lemma\EE For Q = %[F[<f>* + f ) + F((f>* - | )] and v such that 



(2.33) 



cos Z 



1 b 



system (I2.30P has an equilibrium state at $ = |, ^ = (j)*, R = (Z — 9 2 )/uj, where <p* is 
given by fT2TT3l (see also fl2TTl) . fl2TT6l) ) By fl2TT3p . fl2TTl) . the linearization matrix at such 
equilibrium is 

- Pl 

-p 2 cv 

-\e- a ^ z - e ^^smZ + \cv 
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where p x := ae~ a(z_92)/a; (sin Z cos 6 - cos Z sin 6), p 2 := bue' a{z ~ e2)/u {acosZ + usinZ). 
This matrix has three zero eigenvalues at v — v* provided 

(2.34) D{u*) := Pl p 2 + as*{~a>* - h - e -^ z -^)l" s i n Z) = 0. 
At 7 ^ system fl2T34l) . fl2T33|) for Z = Z* transforms into 

cos 2 Z* [l-4 7 2 o^ sin fl]+2 7 sin Z* cos Z* [l+2 7 ^(a cos 6-u sin #)]+4 7 2 ^w cos 6 sin 2 Z* = 0, 

and it is easy to check that the solvability of this equation is given by condition (12.181) . 
Moreover, solutions satisfy 

(2.35) D\v*) ^ 0. 

If 7 = 0, condition (I2.33f) gives cosZ* = 0, and one may check that condition (12. 18ft in 
this case guarantees the solvability of equation (12. 34ft for v* and the fulfillment of ( 12. 35ft . 
It follows from (I2.17P that px j2 ^ at the solutions (hence u* ^ 0) and that 

(2.36) acos(Z* -6) +wsin(Z* -6) ^ 0. 

At v — v* (the triple zero bifurcation moment) the vectors 

/ -cv Pl \ / \ 

Vi = , v 2 = \ cv , v 3 = 

V "PIPS / V / 

form a Jordan base. At v close to v*, take Z satisfying (I2.33f) and close to Z*, and denote 

/ R- (Z-OJ/u \ ( -cu Piyi 

(2.37) $ - f = yxvx + j/ 2 t; 2 + 2/3^3 = cz/ y 2 

V * - 0* / V 2/3 - PlP2 J/1 

System (I2.30p takes the form 

yi = 2/2 + (y 2 ), 

(2.38) y 2 = y 3 + 0(y 2 ), 

m = £i + e 2 y 2 + pyf + 0(M 3 + 1 2/1 1 ( 1 2/2 1 + |y 3 |) + vl + vl), 

where p = i(p lP2 ) 2 z/[F"(0* + f)+F"(0*-f)] ^0, and £l = §*W* + i )+F(0*-f 
e 2 = D(i/), i.e. (ex, e 2 ) are small parameters which are related by a diffeomorphism to the 
original parameters v and Q near the triple zero bifurcation moment (see (I2.34p . (l2.35p ). 
Scale the parameters as follows: 

(2.39) ex = --s 6 , e 2 = Es 2 

P 

for a sufficiently small s, and for some bounded E. By scaling the time and the variables: 

r -»■ cr/s, yi -»■ Fei/s 3 , y 2 ->• ^Vi/s 2 , 2/3 -> i^i/s, 
we bring system ( 12. 38ft to the form 

(2.40) Y'" = 1 - Y 2 + EY' + O(s) 

(where ' denotes the differentiation with respect to the new, slow time a). 



14 



D. TURAEV AND S. ZELIK 



The limit equation 
(2.41) Y'" = 1 - Y 2 + EY' 

has two hyperbolic equilibria: + : Y — 1, with a one-dimensional stable manifold 
and a two-dimensional unstable manifold W", and 0_ : K = —1, with a two- 
dimensional stable manifold and a one- dimensional unstable manifold W™. At E < 3 
these equilibria are saddle-foci, i.e. each of them has a pair of complex characteristic 
exponents. By |KuTs76j . equation ( 12.4 ip has, at E = E* = — -37==, a solution 

Y(t) = - -tanh(^/l 1/120 t) + — tanh 3 (^/ll/120 t) 

which connects the saddle-focus 0_ with + . This solution corresponds to a curve T |_ 

along which the one-dimensional manifolds Wt and coincide. By [IbRo05], at the 
same E there exists another heteroclinic curve, which corresponds to a transverse 

intersection of the two-dimensional manifolds and Wt- By the transversality, the 

heteroclinicnic orbit T H persists for all E close to E*. The other heteroclinic orbit, F h , 

splits as E varies, and this results |Byk99 BykOO] in the sequence of values Ek — > E* which 



correspond to the existence of homoclinic loops to the saddle-foci + and 0_ (equation 
(I2.4ip is time-reversible, so homoclinic loops to the both saddle-foci appear simultane- 
ously). One can view the one-parameter family f !2.4ip as a smooth curve in the space 
of smooth flows in M 3 ; then the parameter values correspond to the intersections of 
this curve with smooth codimension-one surfaces filled by systems with a homoclinic loop 
to, say, the saddle- focus 0_. Importantly, these intersections are transverse. Therefore, 
fixing any arbitrarily large k, we will have at some E close to Ek a homoclinic loop to a 
saddle-focus close to 0_, for any one-parameter family which is sufficiently close to (I2.4ip . 

Thus, given any sufficiently large k, at E = Ek + O(s) equation (12.40P has, for every 
sufficiently small s, a homoclinic loop r* to the saddle-focus 0_ at Y = Y-(k,s) = 
— 1 + 0(s). Denote as £1,2,3 the characteristic exponents at the saddle-focus, £1 > 0, 
Re £ 2 = Re £ 3 < 0, Im £ 2 = - Im £ 3 ^ 0. As & + £ a + £ 3 « here (the limit equation (l2~4Tp 
is volume-preserving), the Shilnikov condition of chaos, £1 + Re£ 2 > 0, is automatically 
fulfilled. Hence, by |Shi65|, IShi70] we obtain an open region in the parameter plane which 
corresponds to a chaotic behavior (i.e. to the sought basic hyperbolic set A) in equation 
(I2.40p and, equivalently, in the original system (I2.30p . 

To finish the proof we need to show that the set A can be chosen in such a way that 
it will contain a pair of periodic orbits for which (I2.3ip is satisfied. According to Remark 
I6.6[ it is enough to check that the integral of the function g — g\ along the homoclinic 
loop to the saddle-focus 0_ is non-zero. In order to verify this condition, let us rewrite 
the function g in the new variables (Y, Y', Y")\ 

g(Y, Y', Y") =v + ^ e -^-e 2 )/u CQ ^ Z _ ^ + c s s y + ^^e^ 

where C := - ^e-^-^^iacosiZ - 6) + tusm(Z - 6)) ^ (see <^M>)- Let Y = 
Y(a\ k, s) be the solution of (I2.40p that corresponds to the homoclinic loop T^; note that 
Y(a; k, s) — > Y_(k, s) exponentially as a — > ±00. Note also that Y(cr\ k, s) = Y(a; k, 0) + 

/+00 
[g(Y(a; k, s),Y'{a; k, s), Y"{a; k, s)) - g{Y^{k, s), 0, 0)] da = 
-00 
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/+oo 
(Y(a;k,0) + I) da + 0(s 4 ). As k — > +00, the homoclinic loops of equation 
■00 

(I2.4ip approach the heteroclinic cycle T + _ U T_ + U 0_ U + at E = E*, so the homoclinic 
loop Tq to 0_ : {Y = — 1} spends at large k a large time in a neighborhood of the other 
equilibrium, + : {Y = +1}. Therefore, the integral of (Y(a; k,0) + 1) tends to +00 as 
k +00. Thus, for sufficiently large k and sufficiently small s, 

/+00 
[g(Y(a; k, s), Y'(a; k, s), Y"(a; k, s)) - g{Y^{k, s), 0, 0)] da + 0. 
-00 

By Remark l6.6[ this proves the lemma, which finishes the proof of the theorem as well. □ 

The proof of the following proposition is standard, cf. [Zel03] . 
Proposition 2.4. Let the non-linearity H satisfy 

(2.43) ReH(z) -z > -C; \H(z)\ < C(l + z 2 ), zeR+ 

for some constant C independent of z. Then for all sufficiently small fi equation (12. ip is 
well-posed in the space C&(R) of uniformly bounded continuous functions and generates a 
dissipative semigroup S(t) t >o in C&(1R) ; and this semigroup possesses a global attractor A. 

The attractor is defined here as follows. Let S(t), t > 0, be a semigroup acting on the 
space Cft(R). A set A C C&(R) is a global (locally-compact) attractor of this semigroup if 

1) A is bounded in C;,(R) and compact in Ci oc (R); 

2) A is strictly invariant: S(t)A = A, t > 0; 

3) as t 00, the set A attracts, in the topology of Q oc (R), the images of all bounded 
subsets B C Cf,(R), i.e. for every neighborhood O of A in the local topology and for every 
bounded B C C 6 (R) there is a time T = T(0, B) such that S(t)B C 0(A) for all t > T. 

Remark 2.5. It is well-known (see e.g. [MirZe08, IZel04j ) that, in contrast to the case of 
bounded domains, the global attractor is usually not compact in C&(R) if the underlying 
domain is unbounded. However, attractor's restrictions to every bounded subdomain 
remain compact. The attraction property itself holds, too, in this local topology only. 

A characteristic property of the global attractor is that it consists of all initial conditions 
which give rise to globally defined solutions. Namely, a function u (x) G C&(R) belongs to 
the attractor if and only if there exists a function u(t,x) G /C such that Uq(x) = u(0,x). 
Note that due to the invariance of the equation with respect to temporal and spatial 
translations, the boundedness and local compactness of the attractor mean also that the 
set K, of the solutions which are defined and bounded for all (t, x) G R 2 is bounded in 
C7 6 (R 2 ) and compact in C ioc (R 2 ). 

Thus, we may define the space-time entropy of the attractor as the space-time entropy 
of the set K: h s _ t (A) := /i s _ 4 (/C) (see (12. 2p ; more discussion and a comparison with other 
definitions can be found e.g. in (MirZe08llZel04j ). As we mentioned (see |CoEc99tlZel07j ). 
the space-time entropy of the attractor of the Ginzburg-Landau equation is finite: 

(2.44) hs- t (A) < 00. 

The next Section gives an explicit example of a scientifically relevant equation with 

(2.45) h a - t {A) > 0. 
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3. Attractor of positive space- time entropy in a perturbed nonlinear 

Shrodinger equation 

Here we prove the following 

Theorem 3.1. Given any sufficiently large (3, there exist (continuously depending on j3) 
intervals of values of S, p, E\ > 0, e 2 and p such that the attractor of the equation 

(3.1) d t u = (1 + i/3)d%u - (1 + i5)u + (i + p)\u\ 2 u - (e 1 + ie 2 )\u\ A u + p 
has strictly positive space-time entropy. 

Proof. The global attractor of equation ( 13. II) exists at E\ > according to Proposition 
12.41 By theorem 12.11 in order to prove (12.451) it is enough to show that the equation 

(3.2) (1 + i(3)d 2 x ll - (1 + iS)U + (i + p)\U\ 2 U - (e 1 + ie 2 )\U\ A U = 

has a non-degenerate symmetric localized solution at some S that depends on the other 
parameters (3, pi j2 , £1,2, and that conditions (I2.14l) . (l2.17p . (l2.18p are satisfied at 7 = 0. The 
localized solution of the ODE (13.21) corresponds to an intersection of the two-dimensional 
stable and unstable manifolds of the hyperbolic equilibrium at U = 0. Because of the 
phase-shift symmetry, when these manifolds intersect they coincide. The soliton non- 
degeneracy conditions imply (among other things) that as 5 changes the manifolds split 
with a non-zero velocity. It follows that a non-degenerate soliton will persist at small 
perturbation of the nonlinearity, provided a small adjustment to the value of 5 is made 
(see more in |AfMi99, AfMiOlJ). Thus, it is enough to consider the cubic equation 

(3.3) (1 + ip)%U - (1 + iS)U + + p)\U\ 2 U = 0; 

once the existence of a non-degenerate soliton is established for this equation, it can 

be carried on to the equation (13.21) for all sufficiently small £x,2j an d since conditions 

fl2TlIj) . fl2TT7j) . fl2TT8|) are open, they will persist as well. 

1 B(l-2w 2 )-3w . l-w 2 + 2wB _ 

Let us choose (3 = — , p = 6 = — ^- — — for some small 

B 1 — 2w z + SwB B(l — w z J — 2w 

B > and w such that B > , 2w 2 ■ Then, if we define 

(3.4) U[x) := d 1 U(xd 2 ) 



where d\ = y j^^+3wB ' ^ 2 = a/1 — w 2 — 2^, we obtain the following equation: 

(3.5) (i + B) [d 2 x U - (1 + iu) 2 U + (1 + iu){2 + iu)\U\ 2 U] = 0. 
It has a localized stationary solution (see e.g. |AfMi99t lAfMiOlj ) 

(3.6) U* = . - . . 

The linearization operator £ is given by 

dp := (i + B) [d 2 x ip - (1 + iufip + 2(1 + iw)(2 + zw)|C/*|V+ 

+ {l+iu){2 + iuj) 2 Ul <p] 

The localized functions 

(3.8) ipi(x) = -d x U*(x), <p 2 (x) = iU*{x) 
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(the odd and, respectively, the even one) belong to the kernel of C. 

/+oo 
(p(x)ip(x)dx, so the conjugate to 
-00 

(I3.7P operator is 

£ty == (i + B) [d 2 x i> - (1 + %ujfi) + 2(1 + iu)(2 + iw)|Z7,|V] + 

( °' 9) +(-i + - iw)(2 - zw)(f/*) 2 ^. 

As £ has two zero modes, one even and one odd, the same holds true for the conjugate 
operator C) . At w = B = the equation for zero eigenfunctions of £t reads as 

(3.io) d 2 j) - i> + ay 2 ^ - 2r 2 ^ = o, 

where we denote 

1 



(3 ' U > rW = ch(z)< 

note that 

(3.12) r"(x) = r - 2r 3 , r'"(x) = (i -6r 2 )r'(x). 

It is easy to see that the odd and even localized solutions of (13.101) are 

(3.13) ^x(x) = iY'(x), ^a(x) = T(x). 

We will look for asymptotic expansions of these solutions at small w and B. By ( I3.9h . 
the localized zero modes of C) satisfy 

ij/'(x) + AT 2 i) - 2T 2 i) = 

( °' ' l) = iw [(2 - QT 2 )iP - (3Y 2 + AT 2 In Y)$] + 4i£r 2 ^ + 0(w 2 + fi 2 ) 

(we take into account that U* depends on w as well: by fl3.6IUI3.lip U 2 = r 2 (l — 2iw\nT + 
0{w 2 )), while |[/,| 2 = T 2 ). By (ETT1) . we have 

(3.15) ip = u + iv + 0(w 2 + B 2 ), 

where 

u"{x) -u + 2T 2 u = -wv{2 - 3r 2 + AY 2 In Y) + ABvY 2 , 
v"(x) -v + 6Y 2 v = wu(2 - 9r 2 - 4Y 2 lnT) + ABuY 2 . 

By fl3.15IUI3.16p . fl3.12p the two sought localized solutions of (13.141) are given by 

(3.17) ^ 1 =tY'(x) + S(x)+0(w 2 + B 2 ) } ip 2 = Y(x) + iQ(x) + 0(w 2 + B 2 ), 
where S and Q are real, decaying to zero, as x — > ±oo, functions which satisfy 

(3.18) S" - S + 2Y 2 S = -w(2 - 3Y 2 + AY 2 In Y)Y'(x) + ABY 2 Y'(x), 

(3.19) Q" - Q + QY 2 Q = w(2Y - 9Y 3 - 4r 3 In Y) + ABY 3 . 

To find S(x), we multiply (I3.18P to Y(x). The equation will take the form (see f)3.12p ): 

YS" - Y"S = -wiY 2 - r 4 + T 4 In r)' + B(Y A )'. 
By integrating this equation with respect to x, we find 

rs" - y's = -w(y 2 - r 4 + r 4 in r) + by 4 



(3.16) 
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(there is no integration constant in the right-hand side, since both S and T tend to zero 
as x —>■ ±00). By solving the first-order equation, we finally obtain 

S(x) = -wT(x)(x - J T 2 dx + J T 2 lnTdx) + BT(x) J T 2 dx 
<3 ' 20) --^(2*-^(2-bch(*))) + B#i 



ch(x) ch(x) 1 ch (x) 

Similarly, by multiplying ( 13.1 9p to T'(x) and integrating the obtained equation, we find, 
with the use of (13TT2|) . that T'Q' - T"Q = w(TT" - T 4 In T) + BY 4 . The solution is 

Q(x) = w(xT'(x) - T - T'(x) [ ^^dx) + BT'(x) [ -^dx 

(3.21) w J ( P ) J (H 2 

= rr— -(2a;sh(x) + ch(x) + ch(a;) lnch(x)) + i?ch(a;). 

ch (x) 

It is immediately seen that functions 5* and Q given by (I3.20I) . (I3.21I) are localized indeed. 
Moreover, S is odd and Q is even, so by plugging (13.201) and (13.211) in (I3.17p . we obtain 
the odd (ipi) and even (tp 2 ) zero eigenfunctions of £' '. 
One can also compute (see (13.81) ) that 

r+00 

Re / tpi{x)tf\{x)dx = 



(3.22) ^ J -™ 

/* + 00 /* + GO 



/-\-OQ /»H~CX3 r\ 

ST'dx + w (r') 2 (l + lnV)dx + 0(w 2 + B 2 ) = -B + 0(w 2 + B 2 ), 
00 J —00 " 



Re / i/j 2 (x)(p2(x)dx = 
(3.23) ^ J -°° M 

\ / /"+oo r+00 

= -w T 2 \nTdx- QTdx + 0(w 2 + B 2 ) = 2(2w-B) + 0(w 2 + B 2 ). 
J —00 J— 00 

As we see, these inner products are non-zero for the values of B and w that we consider 
here (small B,w such that B > and i? > 1 2 ^^ ). This shows that there are no adjoint 
functions to the eigenfunctions (13.81) . The absence (at small B,w) of eigenvalues on the 
imaginary axis follows from |KaSa98j . Thus, the pulse U = U*(x) is non-degenerate. 

Returning to the non-rescaled variables, we find that the soliton U = d^U^x / 'd 2 ) of 
equation ( 13. 3 p is non-degenerate. The corresponding eigenfunctions of C\j are given by 

Mv) = 2B + 0( l ^ + B^ r{x/d2) + S{x/d2) + ° {W2 + 52)) ' 

M*) = Trrrr- J' n( 2j . R2 „ ( r (^ 2 ) + iQ{x/d 2 ) + 0(w 2 + B 2 )) 

2d 2 ({2w — B) + 0(w z + B 2 )) 

(we normalize them so that (EE]) is fulfilled, see (j3T22|) . (l333|) ). By d3~20l) . (I3T2T1) . we find 

. , > 3d 1 (i-w(2-\n2)-B + 0(w 2 + B 2 )) _ (1+iw)M/da . , , 
M*) B + o{w2 + B2) sign(x), 

M ^(1 +iB- iw{l + In 2) + Q( W 2 + B 2 )) (1+tw)lxl/d2 
^ X) ~ d 2 (2w - B + 0(w 2 + 5 2 )) 
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as x —> ±00, so lo = —w/d 2 , a = l/d 2 , and the coefficients s and q in (I2.10p . fl2.lip are 

3tji(j - w(2 - In 2) - B + Q(w 2 + B 2 )) 
S ~ B + 0{w 2 + B 2 ) ' 

_ d x d 2 {l +iB- iw{l + In 2) + 0(w 2 + B 2 )) 
Q ~ 2w - B + 0(w 2 + B 2 ) ' 

It is easy to see that all conditions fl2TT4j) . fl2~T7j) . fl2~18|) hold at small w ^ 0, B > 0. □ 

4. Normally-hyperbolic manifolds for lattice dynamical systems 

In this and the the next Sections we study a class of lattice dynamical systems which 
includes systems describing weak interaction of solitons localized in space and chaotic in 
time, e.g. system fl2.28l) . fl2.29l) . We start with a skew-product system of ODE's 

(4-1) y'(t) = f(y), p'{t)=g{y), 

where /, g are C r , r > 1. We assume that y e R n , p G M m ; for more clarity we will denote 
the space of y variables as Y and the space of p variables as P. We will further assume 
that the y-part of our system: 

(4.2) y' = f(y), 

possesses a bounded, uniformly-hyperbolic invariant set A. 

Recall that the hyperbolicity means that for every point of A there are two subspaces, 
N s (y) and N u (y), such that the following holds: 

1) N s (y) and N u (y) depend continuously on y G A, 

2) the direct sum of N s (y), N u (y) and N c (y) := Span(y) := {\f(y)\\ G M} constitutes 
the whole of R n , 

3) given any orbit y(t) from A, each of the families of subspaces N s (y(t)) and N u (y(t)) is 
invariant with respect to the flow of system (14. 2 p linearized about the orbit y(t), 

4) the linearized flow is exponentially contracting in restriction onto N s (y(t)) as t — > +00 
and in restriction onto N u (y(t)) as t — > —00 (the flow, then, is expanding on N s (y(t)) as 
t — > —00 and on N u (y(t)) as t — >■ +00). 

The linearized system is 

(4-3) jv = f(y(t))v. 

Since v(t) = y(t) = f(y(t)) is a uniformly bounded solution of it, there exists a uniformly 
bounded non-zero solution y*(t) for the conjugate system 

(4-4) j t v = -f(y(t)) T v. 

As y*(t) solves IIP]) , it follows that f t (y*(t) ■«(*)) = for every solution v(t) of (Oil , 
i.e. (y*(t) • f (£)) stays constant. Thus, since the solutions of (14. 3 p which lie in N s (y(t)) © 
N u (y(t)) tend to zero either as t — > +00 or as t — > —00 and y*(t) is bounded, we find that 
this constant is zero for every v G N s (y(t)) © N u (y(t)), i.e. the vector y*(t) is orthogonal 
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to N s (y(t)) © N u (y(t)) for all t. This condition defines y* up to a scalar factor; we fix it 
by normalizing y* in such a way that 

(4.5) (y*(f) • y{t)) = 1. 

The exponential dichotomy for system (14 . 3 [) restricted to v (t) G N s (y(t)) © N u (y(t)) 
implies that the equation 

|«(«) - /'(*(*))« = MO 

has a unique uniformly bounded solution v(t) G N s (y(t)) © N u (y(t)) for any uniformly 
bounded function G N s (y(t)) © N u (y(t)). It is more convenient for us to express this 
property in the following equivalent way: the equation 

(4.6) j f v(t) - f'(y(t))v + (y*(t) ■ v) y(t) = h(t) 

has a unique uniformly bounded solution v(t) given any uniformly bounded function h(t). 
More precisely, equation (14. 6p defines a linear operator L y : h h-» v such that 

(4.7) |M| <C A \\h\\. 

The assumed uniform hyperbolicity of the set A means that the constant C\ in (14.71) can 
be taken the same for all orbits y G A. 

Take a countable set of equations of type (14.11) . This produces an uncoupled LDS 

(lattice dynamical system): 

(4.8) y' k (t) = fk(y k ), p'k(t)=g k (yk), kez 

We assume that the derivatives of fk and gt up to the order r are uniformly continuous 
and bounded for all k, and that for each k the k-th individual ODE's in the LDS has a 
hyperbolic set A/-, all these sets are uniformly bounded and uniformly hyperbolic for all 
k (the uniform hyperbolicity means in our approach that the constant C\ in (14. 7p can be 
taken the same for all k). In the example considered in Section [21 the individual ODE's 
are identical to each other, so the uniformity with respect to k holds trivially. 
By introducing Banach spaces 

Y:=/ 00 (F), ||y|| ¥ := sup keZ \\yk\W, y := {yk}kez, 
P:=/oo(P), ||p||p := sup fceZ ||pfc||p, P := {PkjkeZ, 

we may write the LDS as 

y'(t) = f(y), P / (t) = g(y), 

where f := {f k }kez, S ■= {9k]k&- 

The subject of our study will be a coupled LDS, obtained by a small smooth perturbation 
of this system. Namely, we consider 



(4.9) 




f(y) + eF e (y,p), 
g(y)+eG £ (y, p) 



where e is a small parameter, and F e and G e are C r -functions YxP -> Y and, respectively, 
Y x P — > P; by "C r " we mean, here and below, that all the derivatives up to the order 
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r exist, are uniformly continuous and uniformly bounded. We also assume continuity (in 
C r ) with respect to e, so 

(4.10) ||F e ||or+||G e ||or <C 

where C is independent of e. 

Let y°(t) := {y k (t)} k£Z be a sequence of arbitrary orbits y k {t) G A k ; we will say that 
y°(i) is an orbit from A°°. Each orbit y k (t) defines a curve in the Y-space. The direct 
product of these curves, times the space P, is a C r -submanifold of Y x P, we will denote 
is as Wy - Given an orbit y°, the corresponding manifold W y0 is given by the equation 

(4.11) Vk=Vk(<l>k), kez, 

where the "phases" (f) k run all real values, independently for different k. If we introduce a 
Banach space \1/ of the bounded sequences $ := {<fik\kei with the uniform norm ||$|| := 
sup fcgZ \4>k\, then Wj is a C r -embedding of ^ x P into Y x P. Obviously, W y0 is invariant 
with respect to the non-coupled LDS (14. 8p . Moreover, this manifold is normally- hyperbolic 
(as each of the orbits y k is uniformly-hyperbolic). It is a well-known general principle 
that normally-hyperbolic invariant manifolds persist at small smooth perturbations (see 
|Fen71l IHiPS77j ). The next theorem shows that this principle holds true in our setting. 

Theorem 4.1. For all sufficiently small e, given any orbit y° G A°° there exists a uniquely 
defined C r -manifold W y o e C Y x P, which is invariant with respect to system (j4.9j) . 
depends continuously on e (in C r , uniformly with respect to y°), and coincides with W y0 
at e = 0. Namely, W y o £ is given by 

(4.12) y k = U*($,p,e) := y° k ((/> k ) + V fc ($, p, e), 
where 

(4.13) \\V k \\cr-i = 0(e), \\Y k \\ C r = o(l) e _ , 
uniformly for all k G Z and all y° G A°° . 

Proof. We start with some preliminary constructions. Define the exponential a-norm 
\\h\\ a := sup tgM e~ a '*' \\h(t) \\ on the space of continuous, uniformly bounded functions h; 
e.g. || ■ 1 1 o is just the C°-norm. 

Lemma 4.2. For all small a > and v > 0, for all functions <p(t) such that 

(4.14) \(f)\t) -l\<u for allt Gl, 

and for any A(t) and b(t) sufficiently close (in C°) to f'(y(t)) and, respectively, to y*(t), 
the equation 

(4.15) j t v(t) - A{ct>{t))<t) + (Km ■ v(t)} f(y(m) = Ht) 

is uniquely solvable for any uniformly bounded function h(t), and the corresponding linear 
operator L^-.hv-^v satisfies 

(4.16) \\v\\ a <C, 
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Moreover, the operator is Lipshitz with respect to <fi: ifvi(t) andv2(t) are the solutions 
of equation f)4. 15j) which correspond to two different functions <fii(t) and 02 (t) (and to the 
same right-hand side h), then 

(4.17) \\V 2 - Ui||a < K||/i|| o ||02 - <Pl\\a 

for some constant K, proportional to the C x -norms of A, b and f. 

Proof. A uniformly small continuous perturbation of the time-dependent coefficients in 
the left-hand side of (14.61) does not destroy its unique solvability property. Hence, equation 

(4.18) jv{t) - A(t)v + (b(t) ■ v) f(y(t)) = h(t) 

has a unique uniformly bounded solution v{t) given any uniformly bounded function h(t); 
moreover, for the corresponding operator L : /i i— )■ v estimate (14 .7p holds (we assume that 
the constant C\ in (14.71) was taken with a margin of safety, so all our small perturbations 
of the equation do not change C\). Note also, that given any function <p(t) that satisfies 
(I4.14p . if we introduce a new time r = 0(t) in the equation (14.151) and a new function 
Vnew by the rule v new ((fi(t)) = v(t), then the left-hand side of equation will be 0(z/)-close 
to the left-hand side of (14.181) . For sufficiently small v this gives us the unique solvability 
of (I4.15P and estimate (14.161) with a = 0. 

Next, we note that a multiplication of the functions v and h in (14.151) to any smooth 
function of t with uniformly small derivative just results in a uniformly small correction to 
A(4>(t)). This immediately shows the unique solvability of equation (14.151) in any weighted 
space with a sufficiently slowly growing weight; e.g. we obtain (14.161) for all small a. 

In order to show the Lipshitz property of with respect to 0, we note that 

v 2 -v x = {(Afa) - A^))^ - [(&(&) • Vl ) f(y(<f> 2 )) - (6(0i) • «i> fWi))]} ■ 

Now, since A{4>), b(<p), f(y(4>)) are smooth - hence, Lipshitz - with respect to 0, and since 
v±(t) is uniformly bounded by (14.71) . we immediately get (I4.17P from (I4.16p . □ 

Further we will use 

/+oo p+oo 

y*{t + sfi)as)ds, A{t) = / f(y(t + sfi))Z(s)ds, 

-oo J — oo 

where /i is a small constant and £ > is such that f_°° £(s)ds = 1. At fi = we have 
b = y* and A = f'(y); at small /i the functions A{t) and b(t) are close, respectively, to 
f'(y(t)) in C r ~ l and to y*(t) in C r (we have y*(t) e C r as it satisfies equation (I4.4p ). 
Thus, uniformly for all t, we have 

(b(t).f(y(t)))-l:=c(t) = 0(»), 

(4.20) 

b'(t) + f'(y(t)) T b(t) = o(l)^ , A(t) - f'(y(t)) = o(l) r+o 

(see (@3D,(|42D). By taking ^ G C°° and such that f+™ \£{s)\ds < oo, we will make A(t) 
and b{t) at /i ^ more smooth then f'(y(t)) and, respectively, y*(t), namely we will use 
A which is at least C r and b which is at least C r+1 ; the price is that the last derivatives 
do not stay bounded as fi — > 0, however we have en estimate 

(4.21) Holler = Ofr- 1 ), \\b(t)\\ C r + l = OQT 1 ). 
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The next proposition describes the way we coordinatize a small neighborhood of the 
curve w% : y — y k (t) in the F-space (y k (t) is an orbit from the hyperbolic set A k of the 
k-th subsystem of the uncoupled LDS ( 14. 8p ). Let w k : y = y k (t) De a curve, 7-close to w k 
on some, finite or infinite, interval I of t, i.e. there exists a smooth time-reparametrization 
i/i(t) such that \\y° k (il>(t)) - y k {t)\\ Y < 7 at t G I. 

Lemma 4.3. There exists 7 > (independent of the choice of the orbit y k G A k and 
independent of k) such that if j < 7, then there exists a uniquely defined on I function 
4>(t) such that 4> = ip + 0(7), ^ = 1 + o(l) 7 _ > o? an d 



where b k is given by (14.191) at some small \i. 

Proof. The derivative of the left-hand side of ( I4.22p with respect to at constant y k is 



(b'M ■ (y k - vim) - H<P) ■ y°M) = 0(y k - y°M) - {h{<p) ■ f k (y° k (4>))) ■ % < SM, 



it is bounded away from zero, provided y k — yl((fi) is sufficiently small. Thus, by the 
implicit function theorem, for any point y k from the (sufficiently small) 7-neighborhood 
of the point y^.(<p), we have a uniquely defined 4>{yk) which satisfies (I4.22p and condition 
4>{y^{ip)) = ip- Moreover, depends smoothly on y and the derivatives are uniformly 
bounded. So, as — j/&(i)||y < 7, we also have \\ip(t) — 0(t)||y = 0(7), as required 

(we denote <f>(t) := (j>(y k (t))). □ 

Let us now proceed to the proof of the theorem. Let y°(i) be an orbit from A°°. The 
sought invariant manifold W y o e consists of all solutions of the LDS (14. 9 j) which stay for all 
times in a small neighborhood of the manifold W° . This means that, for every fceZ, the 
k-th component of y(t) stays uniformly close to the corresponding curve w k ) :y = y k ((fi) 
in the F-space. In other words every trajectory (y(t),p(t)) G W y o j£ satisfies 



where the functions v k (t) are uniformly small. By Lemma |4.3[ we may always assume 
that the parametrization is chosen so that (14.221) is fulfilled. By differentiating (14.221) 
with respect to t we get 



where we denote z k {t) := y^((p k (t)). By multiplying both sides of this equation to b(cj)k(t)), 
and taking (I4.22p . (l4.24l) and (I4.20p into account, we obtain the following equation for the 
evolution of cf) k \ 



(4.22) 




(4.23) 



y k (t)=y° k (<f> k (t))+v k {t), kez, 



(4-24) (h(Mt)) ■ v' k (t)) = -<j>' k {t) (b' k (Mt)) ■ v k (t)) . 



Now, plugging (I4.23P into the first equation of (14.91) gives 
(4.25) v' k {t) + <t>' k {t)f k {z k {t)) = f k {z k (t) + v k {t)) + eF £ik {z + v(t), p(t)), 



(4.26) 



(j)' k (t) = l + g fe (v,$,p), 



where 



(4.27) q k : 



([fk(z k + V k ) - /fc(-Zfc)] • b k ((j) k )) + (&fc(0fc) • Vk) + £ (FeA z + v > P) • ftfc(</>fc)) 

1 + c k {(j) k ) - (b' k ((f) k ),v k } 
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Equation for the v -components can now be obtained by plugging (14.261) into (14.251) : 

(4.28) v' k {t) - A k (<f) k )v k + (b k {<f) k ) ■ v k ) f k {z k ) = Q k (y, $, p) - q k (v, $, p)f k {z k ), 
where 

(4.29) Q k := f k (z k + v k ) - f k (z k ) - A k ((p k )v k + eF £tk (z + v, p). 
Equation for the evolution of p(i) is given by the second equation of (14.91) : 

(4.30) p' k (t) = g k (z k + v k ) + eG £ik (z + v, p). 

We remark that if we choose /J = in (14.191) . equations (14.271) and (14.291) are simplified 
and reduce to 

(4.31) q k = — ]%T^4rr \ \ ' Qk = fk(z k + v k )-f k (z k )-f k (z k )v k + eF £jk (z + v,p) 

1 + {ykWk),f k {z k )v k ) 

(see (I4.20p ). However, the functions q k and Q k will be only C r ~ l with respect to $, 
therefore we use small non-zero /x - in order not to lose the last derivative (and to be able 
to treat the case r = 1). 

By multiplying both sides of P~2"g|) to b k (<j) k (t)) and using (jOEl) . (H~20J) we find that 

j MMt)) ■ v k (t)) + (h(Mt)) ■ v k (t)) (1 + 0(ji)) = 0. 

This equation has only one bounded solution: (b k (<p k (t)) ■ v k (t)) = 0; therefore, since b 
is uniformly bounded, we find that every uniformly bounded solution v(t) of the system 
(14. 26|) . (14. 281) . (14. 301) (with k running all integer values) satisfies (I4.22j) . Hence, it satisfies 
(I4.25p . Thus, the solutions of system (I4.26[) . (14. 281) . (14. 301) whose v(t)-component is uni- 
formly small give us all the solutions of system (14.91) which stay uniformly close to the 
manifold W° (i.e. all the solutions which comprise the sought invariant manifold W y o e ). 
We show below that for all small 6 > the solution of (TC26]), (j428]), (TOO]) for which 



(4.32) IK(*)l|y < 8 (keZ, teR) 

exists and is defined uniquely for any given initial condition $(0) and p(0). 

In order to prove the existence and uniqueness of the (small v) solution, we will show 
that it can be obtained as a fixed point of a contracting operator on an appropriate 
space. Namely, we consider the set Xg,u of all functions (v(i), <&(£), p(t)) belonging to 
space C loc (R,Y x * x P) such that (j021> and (fl~T4"j) hold for all k and t for which the 
following norm is finite: 

(4.33) ||v,$,p|U= sup e- a ^m a x{\\v k (t)\\,\Mt)\,4Pk(t)\\}, 

fcez,ieR 

where a > 0, and k > is assumed to be sufficiently small. Obviously, the set X^ v is a 
complete metric space with respect to that norm. 

Note that in the limit lim £ _;.o,v-^o the functions q k , Q k given by (I4.27p . (l4.29p tend 
uniformly to zero for all fceZ, and in the limit lim^+o lim e ^o,v^o their first derivatives 
with respect to v, $ and p tend uniformly to zero too (see (14.101) . (I4.20p . (l4.2ip : the order 
of the limits is important: e and v first, then //). The first derivative of the right-hand 
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side of (I4.30p with respect to p is also uniformly small. Thus, if we rewrite system 



where ($°, p°) G P x $ is arbitrary and the operator L is defined by equation (14.151) . then 
it is easy to check that the right-hand side of (I4.34[) (for every fixed $° and p°) defines a 
contracting operator T : X$,v ~^ X$ )V for every exponential norm with a sufficiently small 
weight aQ. Namely, we first fix small v and such that the operators L^ k will all be 
defined and Lipshitz with respect to 4> (see (I4.17p . actually, we may fix v of order e); the 
operator of integration J^(-)dt is also Lipshitz in the ao-norm, with the Lipshitz constant 
— ; then we choose a sufficiently small fi for which the Lipshitz constants of qu and Qk 
can become small enough as e and v tend to zero; then we see that one may choose k 
sufficiently small such that for all sufficiently small e and 8 the Lipshitz constant of the 
right-hand side of (14.341) on the space Xg jU is less than 1, which means the operator T 
is contracting indeed (we need to introduce the small factor k in the definition of norm 
on Xg tU because the derivative of gk with respect to Vk and 4>k, though bounded, is not 
necessarily small). As at e = and v = the v-component of the image by T vanishes, 
the contractivity of T implies that given any small S the condition (14.321) is invariant with 
respect to T for all sufficiently small e; i.e., the TXs,u C Xg v . 

By the Banach principle, the iterations by T of any initial element from Xg )U converge 
to a uniquely defined limit in Xg }U , the fixed point of T. Thus, we have shown that every 
solution which stays sufficiently close to the manifold W y0 for all times can be found as the 
uniquely defined solution of (I4.34j) . Therefore, the union of all such solutions comprises 
the sought invariant manifold W y o £ given by (14.121) where the function is defined 
by the map ($(0),p(0)) ffc(O). Note that this map (hence the manifold W y o £ ) is 
Lipshitz continuous, since the contracting operator T is Lipshitz continuous with respect 
to ($°, p°). We omit the proof of the smoothness of this map, as it is completely standard 
(yet laborious): one may show that the operator T is smooth on a scale of Banach 
spaces corresponding to different weighted a-norms on Xg iV (cf. [Gi Va87[ IShSTCOT] ). or 
alternatively check, by fiber-contraction arguments, that the iterations of by T of an initial 
element from Xg }U converge to the fixed point of T uniformly along with the derivatives 
with respect to ($°,p°) (cf. [Mie86l [Mie58] ). 

In order to finish the proof of the theorem, it remains to show estimate (I4.13p . The 
C'-part is obvious, as v = solves (14.341) at e = 0, and the fixed point of a contracting 
operator which depends on a parameter continuously must depend on the same parameter 
continuously. To show the C"" _1 -estimate, we note that when the right-hand side of (14.341) 
depends smoothly on some parameter, the solution must also be smooth with respect to 



mB,<SM,<SM as 




(4.34) 
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the same parameter. In particular, if we rewrite system ( 14. 9 p as 

y'(t) = f(y) + aF £ (y, p), p'(f) = g(y) + aG £ (y, p), 

then v(i) will depend C r -smoothly on a as well, which immediately gives (I4.13P if we note 
that v(i) = at a = and plug a = e back. □ 

Remark 4.4. The theorem also remains true if the number of systems coupled in the 
LDS is finite, i.e. if the index k runs a finite set instead of Z. Then the range of e values 
for which the corresponding invariant manifolds exist will be independent on the number 
N of systems in the LDS - provided the constant C in the bound (14. lOj) on the norm of 
the coupling terms is independent of N. Note that condition (14. 10p does not requires that 
the coupling is local, it just means that the "total coupling strength" for each subsystem 
in the LDS is bounded independently of the total number N of subsystems involved. 

Remark 4.5. Given any symmetry in system (14. 9p . if the set A°° obeys the same sym- 
metry, then system of invariant manifolds W y o e inherits the symmetry for all small e - by 
uniqueness. A basic example of such symmetry is invariance with respect to spatial trans- 
lation k — > k + 1 (in this case, the coupling terms F, G are shift-invariant, the individual 
ODE's (14. 8 p are the same for all k, and the sets should be chosen the same). 

Theorem 14.11 allows us to construct a huge number of special solutions of the weakly 
coupled LDS ( 14. 9p . Indeed, for every y° G A°°, one can construct the associated manifold 
W y o >£ and then, for every ($ , p ) G$xP, there exists a solution (y(t), p(i)) in the form 

(4.35) y(t)=v°($(t))+V y o($(t),p(t)) 

where the functions ($(£), p(t)) solve the reduced problem on the center manifold (see 

|V(i) = l + q(V y o($,p),$,p), 

(4.36) I p' = g( y °(0) + V y o($,p)) + £ G £ (y ($) + V y o($,p),p), 

[$(0) = $ , P(0) = p - 

It is interesting to have an expansion in powers of e for the system on the invariant 
manifold. In order to do this we need a sufficient smoothness of the right-hand sides: for 
instance, to find the first order in e approximation to (I4.36P we assume the original system 
to be at least C 2 with respect to all variables and e. In this case we may take fi = in 
formulas (I4.19p . so the function q will be given by (14.311) . As V = 0(e) by (14.131) . we 
immediately obtain the first-order in e approximation to the ^-equation: 

(4-37) <U k (t) = l + s <F M (y°($), p) ■ yUfa)) , 

where yl(s) is the uniquely defined bounded solution of 

^yl(s) = -fM(s)) T yl(s), (yt(s) ■ h(y° k (s))) = l. 

Formula ( I4.37P describes the evolution of phases on the invariant manifolds W and can 
be useful in the study of phase synchronization in coupled chaotic systems (see e.g. 
|RoPK96j l 
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To obtain the approximate p-equation, we need the first-order approximation to v. By 
expanding the first equation in fl4.34j) in e, we find that 

v k = eu k (0) + o(e), 

where the function Uk{t) is given by 

u k = L Mt) [F Jk ,o(y°($(t)),p(f)) - <F fci0 (y ($(t)),p(t)) ■ yl{<t> k {t))) fMiM)))] ■ 

Since, by (I4.36p . (f>k(t) is for all k uniformly 0(e:)-close to (p k (0) + t in the exponential 
a-norm with a > 0, and Pk(t) is for all k uniformly 0(e)-close to p k (0) + f gk(y k ( s + 
(pk(0)))ds, also in the exponential a-norm, it follows from the Lipshitz property of the 
operator (see fl4TT7j) . f l4TT6|) ^ that 

v k = ew k ($,p) + o(e), 

where, given any constant $ and p, we denote as w k (§,p) the value at t = of the 
uniquely defined bounded solution w(t) of the equation 

f t w(t) - f k (y° k (t + <p k ))w + (y*{t + <j> k ) ■ v) f k {y k {t + <f> k )) = 
= J= k - (F k ■ y* k (t + 4>h)) fMtf + fa))), where 
:= F kfl (y°(t + $), p + jl g(y°( S + $))ds). 

By plugging the above formula for v k into the p-equation of (I4.36[) and dropping all 
o(e)-terms, we find that the first-order approximation to the p-equation is 

(4.38) p' k (t) = g k (y M + ew k (<5>, p)) + eG k , (y°(<5> ), p). 

Remark 4.6. Let the uniformly-hyperbolic sets be compact and locally-maximal, i.e. 
there exists 7 > (independent of k) such that for each k every orbit of (14.8j) . which 
stays in the 7°-neighborhood of for all t, belongs to A k itself. Then, for e sufficiently 
small, every solution of the coupled LDS (14.9p whose y-component stays for all times in 
a small neighborhood of A°° belongs to one of the manifolds W y o e . Indeed, given any k 
the fc-th component y k (t) of such solution must be close, after some reparametrization of 
time, to a j-orbit y(ip(t)), which is a countable union of consecutive pieces y{(p)\^[<t>j,4>j+i) 
of orbits from the set A& such that ||y(v?j) — f/i^j ~ 0) II ^ f° r some small 7. It is known 
that when A^ is locally-maximal, any 7-orbit is shadowed by a true orbit, i.e. there exists 
an orbit in A& which is 0(7)-close to y(ip) (after a reparametrization of time). Thus, our 
solution y(t) of the coupled LDS stays for all times close to a (time-reparametrized) orbit 
y° G A°°, i.e. we can write it in the form (|4.23|) . and we showed in Theorem 14.11 that 
every such solution belongs to the invariant manifold W y o e . 

5. A THEOREM ON ASYMPTOTIC PHASE 

In this Section we compare the behavior of orbits of the LDS (14.91) which belong to 
different invariant manifolds W y o >£ . We start with the analysis of the dependence of the 
invariant manifold W y o e on the choice of the trajectory y° G A°°. Clearly, W y o e depends 
on y° continuously: namely, the function V y o e ($,p) (hence - the function U y o e ($,p)) 
in (14.121) is found via an application of the contraction mapping principle, and the corre- 
sponding contracting operator (the operator T defined by the right-hand side of (14.341) ) 
depends continuously on y° in some exponential weighted norm, so on any bounded set 
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of values of $ and p the functions U y o,i )S and U y o,2 )£ will be uniformly close provided the 
trajectories y ' 1 and y 0,2 are sufficiently close in the weighted norm. We need, however, 
a different statement about the closeness of U y o,i e and U y o,2 £ . Note that though the 
manifold W y o e is defined uniquely (as the set of all solutions that for all t stay uniformly 
close to the manifold W y0 defined by ( 14. lip ), the function U y o e in ( I4.12p is defined up to 
an (arbitrary) reparametrization of the space \l/ of phases <j)k- Therefore, when comparing 
functions U corresponding to two different trajectories y° (as we do it below), we should 
describe how the corresponding parametrization choices agree with each other. 

In order to do so we recall the construction used in the proof of Theorem 14.11 Take 
any orbit {y = y°(t)} G A°°; its £;-th component y^it) defines a smooth curve in the 
space Y. Take any other curve w : {y — yk(t)} in Y. By Lemma FO| there exists 7 > 
(independent of the choice of the curves) such that if w stays in the 7-neighbourhood 
of w® for a certain interval of time, then, for every t from this interval, condition ( I4.22p 
defines the projection y1((pk(t)) of the point yuif) G w onto the curve uniquely. We 
will call <pk{t) the phase relative to y°. If we have two orbits, y ' 1 and y ' 2 , from A°°, and 
these orbits are 7-close (7 < 7) on some time interval, then for every curve y = yk{t) 
which stays at the distance less than 7 from both w^ 1 and w®' 2 on this time interval we 
have two phases, 4>\{t) and 4>l(t), relative to y 0,1 and y ' 2 respectively. By Lemma [O] 
(with ip standing for <f>\ and (f> for these two phases are related by a close to identity 
diffeomorphism: 4>l(t) = rjk(4>l(t)) where rj' k ((fi) = 1 + o(l) 7 _>o- For a solution in the 
invariant manifold W y o e the canonical phases (f>k (which we used before) are the phases 
relative to y°. However, if two orbits y 0,1 and y 0,2 are 7-close (7 < 7) 011 some time 
interval, then for solutions in, say, W y o,j j£ the phase ifk relative to y ' 1 is also defined on 
this interval, along with the canonical phase cj>k relative to y 0,2 . 



Lemma 5.1. Let the assumptions of Theorem 4-1 hold. Then there exists 7 > 0, a > 
and C > such that, for all small e, given any T > T > 0, if any two orbits y 0,1 and 
y0,2 j rom satisfy 

(5.1) sup ||y ' 1 (t)-y°' 2 (t)|| ¥ <7, 

te[-T,T] 

for some 7 < 7; then there exists a uniformly close to identity diffeomorphism 77 such that 

(5.2) UIVm^Po) -U y o, 2 (7/($ ),Po)lk < Ce-< T ~ T ^ 
for all p G P and $0 such that 

(5.3) ||$o|k<T . 

Proof. Let (y 1 (t) , p 1 (t)) and (y 2 (t), p 2 (t)) be the orbits on the invariant manifolds, re- 
spectively, W y o,i and W y o,2 such that (y 1 (t) , p 1 (t)) corresponds to the initial condi- 
tion $(0) = $ an d p(0) = p , and (y 2 (t), p 2 (t)) corresponds to the initial condition 
$(0) = 77(^0) an d p(0) = p , where rjk : (fk ^ <Pk is the close to identity diffeomorphism 
which sends the phases relative to y ' 1 to the phases relative to y 0,2 . 

Let 4>k{t) be canonical phases of y l {t) and let v(t) := y x (t) — y°' 1 ( < I ) (t)), so 

(5.4) v(0) = U y o, 1 ($ ,p o )-y°' 1 ($ ) 
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(see (I4.35P ). Let ip k (t) denote the phase of yl(t) relative to y 0,1 , and let u(t) := {u k (t)} k£ - A 
where Uk{t) = yl{t) — Vk {fkit))- By construction, ip k (0) = 0^(0) for all k, so 



a 



(5.5) u(0) = U y o, 2 (7 ? ($ ),p )-y°' 1 ($ 

Denote x\(t) := (v k (t),(J) k (t),pl(t)) and x k (t) := (u k (t), (f k (t),pl(t)); as we showed in 
the proof of Theorem (14. ip the functions x J k (t) satisfy the same system (14.261) . (14.281) . (14.301) : 
x\(t) satisfies this system for all t, while x\(t) satisfies it for all t for which u(t) remains 
small. As the orbit y 2 {t) belongs to the invariant manifold W y o,2 £ , it stays close to 
y 2 '°($ 2 (t)) for all times (where $ 2 is the canonical phase of y 2 ), so by (15.11) the distance 
||u(t)|| between y 2 (t) and its projection to W y0 ,i will remain small for all times such that 
||$ 2 (t)|| < T. Since the time derivative of $ is bounded (see (14.351) ). we have from (l5.3l) 
the required smallness of u(t) for all \t\ < S where 

(5.6) S = O(T-T + l). 

Outside this time-interval we cannot guarantee that the phases <p k (t) are well-defined, 
therefore we modify x\(t) at |£| > 5 — 1. Namely, we consider the functions x\(t) = 
(v k (t) , 4> k (t) , p\{t)) defined by the following rule: 

v 3 k (t) = e (t))u k (t), P l(t) = e (t) P l(t), 



(5.7) 



(t) = e (t)<f k (t) + 9_(t)[M-s + 1) + <f' k (-s + i)(t+s- 1)]- 



+ 9 + (t)[Ms - 1) + <p' k (s - m -S + 1)], 

where 9±(t) are smooth functions IR 1 — > [0,1] such that 6*_(t) equals to 1 at t < —S 
and to at t > —S + 1, while 9 + {t) equals to 1 at t > S and to at t < 5 — 1, and 
9q := 1 — 9 + — 9-. Note that it follows from (I5.7P that f|(t) is uniformly small for all 
t € R since u(t) is uniformly small for all \t\ < S. 

Since both x x (t) and x 2 (t) satisfy system (jOSjl . (Qg|) . flCTD at t G [t x + S,t 2 - S], 
the function x 3 (t) satisfies the same system with a uniformly bounded correction to the 
right-hand sides which is localized at \t\ G [S — 1, S] (and which is denoted below as p). 
Since the initial conditions in $ and p coincide for x 1 and x 2 by construction (recall that 
we choose y x (t) and y 2 (t) such that ip k (0) = 4> k (0)), we find that x 3 (t)i e [_ 00j00 ] satisfy the 
following equation (a perturbation of (I4.34p ) 

; - V fc [Q fc (v 3 , $ 3 , p 3 ) - g fe (v 3 , ^, p 3 )f k (zl) + p kl ] , 



(5.* 



v 



k 



M+t+ [ [g fc (v 3 ,$ 3 ,p 3 )+ Pfc2 ]cit, 
J o 



P 3 k=Pk{0)+ [ [g k {zl + vl) + Pkz] dt + e [ G £ife (z 3 + v 3 ,p 3 )rft, 
Jo Jo 



(k G Z) 



where z\ := 1 {4>l) and the perturbations p k (t) satisfy ||p|| Q0 = 0(e a ° s ). 

Recall that x\{t) = (u k {t),ip k (t),p 2 k {t)) satisfies system (Q5]l . . (ODll at \t\ < S, 
and u(t) is uniformly small on this interval (provided 7 and e are small enough). The 
smallness of u and e implies the smallness of the functions q k in the right-hand side of the 
equation (14.261) for the phases (p k , therefore sup| t | <s \(f' k (t) — 1| is uniformly small for all 
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k. By (15.71) . we find then that sup 4gR \<Pl'(t) — 1| is also uniformly small. This guarantees 
that the operator is defined and Lipshitz in the Q; - norm ( see comments after (14.341) 
in the proof of Theorem 14. ip . Since operator is Lipshitz in the exponential ao-norm, 
and so is the operator of integration J Q (-)dt, we may rewrite (15. 8 p as 

x 3 = Tx 3 + 0(e- aoS ) ao , 

where T is the operator defined by the right-hand side of ( 14. 34 p . i.e. x 1 = Tx 1 , and we 
immediately get that 

Hx 3 -^!!^ =0(e- aoS ), 

since the operator T is contracting (in the norm given by (|4.33p ; note that, as we have 
shown in the proof of Theorem (14.11) . in order to have contraction, both v 1 and v 3 must 
be uniformly small, i.e. must satisfy (I4.32p with a sufficiently small 5, and this property 
indeed holds true when e and 7 are sufficiently small). 

In particular (since v 3 (0) = u(0)), we have ||v(0) - u(0)|| Y = 0(e~ a ° s ), and the lemma 
follows from fl£3} , (ESJ) , . □ 

Remark 5.2. By a shift of time, we obtain that if 

(5.9) sup ||y D - 1 (t)-y w (t)||Y<7, 

tE[Ti,T 2 ] 

then 

(5.10) ||HJ y o,i(^o,Po) -U y °Av(®o), P )||y < Ce- aT ^ 

for all p G P and $ such that for all k the components <fik of $0 satisfy 

(5.11) T x + T <<p k <T 2 -T. 

We may now prove the following theorem, crucial for the next Section. 



Theorem 5.3. Let the assumptions of Theorem \4-l\ hold. Then there exists a > and 
7 > such that for all sufficiently small e and all 7 G (0,7), for every two trajectories 
y ' 1 and y 0,2 from A 00 satisfying the condition 

(5-12) sup||y ' 1 (t)-y°' 2 (t)|| ¥ <7, 

t>t 

given any solution (y 1 (t), p x (t)) from the invariant manifold W y o,i j£; there exists a unique 
solution (y 2 (t), p 2 (t)) from the invariant manifold W y o,2 e such that 

(5.13) Hy 1 ^) - y 2 (i)|| ¥ + Hp 1 ^) - P 2 (t)|| P < C^e^-^, t > t . 
The factor (7(7) tends to zero as 7 — > 0. 

Remark 5.4. Absolutely analogously, for every two trajectories y ' 1 and y 0,2 from A°° 
satisfying the condition 

(5.14) sup||y ' 1 (t)- y °' 2 (t)|| ¥ <7, 

t<t 

given any solution (y 1 (t) , p 1 (t)) from the manifold W y o,i e , there exists a unique solution 
(y 2 {t ) , P 2 (t ) ) from W y o,2 )6 such that 

(5.15) Hy 1 ^) - y 2 (t)|| ¥ + Hp 1 ^) - P 2 (t)|| P < C( 7 )e-^, t < t . 
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Proof. As we explained in Lemma I5.1[ condition (I5.12p (which is an analogue of condition 
( 15.1 p for the case of infinite time interval) implies that for all sufficiently small e, for any 
two solutions (y 1 (t), p 1 (t)) and (y 2 (t), p 2 (t)) from the invariant manifolds, respectively, 
W y o,i ]£ and W y o,2 £ , for all t > to and every k G Z we have well-defined projections of the 
points y\(t) and yl(t) onto the curve y = y^ifk) m the space Y. The position of the 
projection point is defined by its phase ip k , so we have two phases (relative to the same 
orbit y ' 1 ) defined for all t > to'. <p k (t) for the point y\{t) and fl(t) for y k (t). Thus, 

(5.16) yi(t)=yi 1 (fi(t)) + v i k (t), where (b k {<p\(t)) ■ = t > t ; 

here b k is given by ( 14.191) with y* = y^ . 

As the solution (y 1 (t), p 1 ^)) belongs to the invariant manifold W y o,i j£ associated with 
the orbit y ' 1 relative to which the phase is defined, the phases tp\ are just the canon- 
ical phases (f)\. For the solution (y 2 (t), p 2 (t)), as we explained in the introduction to 
Lemma 15.11 the phases ip\ are related to the canonical phases 4>\ by a close to identity 
diffeomorphism r\k : ip\ \-> (fil at t > t ; so, (p\{t) — ^\(t) is uniformly small for all t >t . 

Formula (I5.16P is identical to (I4.23p . (l4.22p . hence (see the proof of Theorem I4.ip 
the functions (v*(t), ^(t), p l (t)), both for i — 1 and i = 2, solve the same system 
( I4.26p .( l4~28|) . (l4.30p (where one should replace <p with (p and z k with y^ l (^k)), for all t 
for which v*(t) remains small. As the solution (y 1 (t), p x (t)) belongs to the invariant 
manifold W y o,i j£ , we have that y l {t) stays close to y 1 '°(<l )1 (t)) for all times, which guar- 
antees the smallness of v 1 . The smallness of v 2 (t) at all t > to follows from the fact 
that (y 2 (t), p 2 (t)) belongs to the invariant manifold W y o,2 e , hence yl(t) stays close to 
Vk (0fc(^)) f° r a ^ times, and because of the uniform closeness of 4>l(t) to (flit) and y k ' to 
2//1' at t > to we obtain the uniform closeness of y\(t) to yl^i^Pkit)) at t >t . 

Thus, we have 

(5.17) Lp\ = 1 + q k (y\ <p\ p J ), j t pl = h k (y\ <p\ p l ) (keZ), 

where qk is given by (I4.27P (the only important thing for us is that qk is uniformly small 
along with its first derivatives), and 

(5.18) h k (v, if, p) := </ fc (yf x (v7fc) + v k ) + eG £jk (y ' l (ip) + v, p). 

By Theorem (14.11) . since the solutions (y*(t), p*(t)) belong to the respective invariant 
manifolds W y o,t )E , we may put 

(5.19) v J (t) = V l (^(t),p i (t)), 

in equations ( 15.171) . where V 1 are certain functions of (<p,p) with the Lipshitz constant 
uniformly small. Namely, the function V 1 is just the function V 1 that defines the manifold 
W y o,i, e by (147121) . while the function V 2 is given by Y(<p, p) = Y 2 (r](ip), p) + y°' 2 (^(v?)) - 
y ' 1 (ip), where rj is the diffeomorphism which sends the phase relative to y ' 1 to the phase 
relative to y 0,2 ; the required Lipshitz property of V 2 follows from the Lipshitz property 
of V and rj. Note that by (H7T2]) 

(5.20) v V P ) = u V P ) - y * 1 ^), vV P ) = u 2 (v(<p), p) - y 0,1 ^), 
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hence, by Remark I5.2[ when (p^ — > +00 uniformly for all k, we have 
(5.21) V%, p) - V 1 ^, p) = 0{e- a '^) 

for some a' > 0. 

It follows from fl5.2ip . fl5.19p . fl5.16p that we will have the required exponential decay 
of || y 1 ^) — y 2 (£)|| if the difference between the corresponding solutions (<p 2 (t), p 2 {t)) and 
(^(t), p 1 ^)) of fl5T7l) .f l5TT9ll tends exponentially to zero as t -> +00. Given (ip 1 (t) , p 1 (t)) , 
the sought, tending to it solution (<p 2 (t) , p 2 (t)) corresponds to the fixed point of the 
operator (</>(£), p(£)) t > (->■ (0(t), p(t))t>o defined by the following equation: 



(5.22) 



/oo 
q(V V, p 1 ), <p\ p 1 ) - q(V 2 (^ + 0, p 1 + p), V l + 0, p 1 + p) 

/oo 
h(VV, P 1 ), </?\ P 1 ) - h(VV + 0, p 1 + p), ip 1 + 0, p 1 + p) 



dt, 
dt, 



where we denote <p{t) := (p 2 (t) — (p 1 (t) , p(t) := p 2 (t)—p 1 (t). Thus, it remains to prove the 
existence and uniqueness of the fixed point of this operator in the space of exponentially 
decreasing functions, and also to show that this fixed point tends to zero as 7 — > 0. 

In order to do this, we first note that because ^(t) and (p 2 {t) grow within linear bounds 
with time, estimate f)5.2ip along with the boundedness of the Lipshitz constants of the 
functions q and h implies that for some a > the operator f 1 5 . 2 2 j) takes exponentially 
decreasing functions (<f>(t),p(t)) = 0(e~ at ) into functions (0(t),p(i)) which are exponen- 
tially decreasing as well, with the same exponent a. 

Recall also that the Lipshitz constant of q is uniformly small (and tends to zero as 
e — > and 7 — > 0). The Lipshitz constant of V with respect to (p is uniformly bounded 
and the Lipshitz constant with respect to p is uniformly small (and tends to zero as e — > 0; 
see fl5.20p . p.l2p .f l4.13|) ). The Lipshitz constant of h with respect to p is of order e and 
the Lipshitz constants with respect to <p and v are bounded (see fl5.18p .f l4.10p ). Thus, for 
the functions under the integrals in f)5.22p . the Lipshitz constants with respect to both 
ip and p in the first equation of f)5.22p and with respect to p in the second equation are 
uniformly small, while the Lipshitz constant with respect to <p> in the second equation 
is uniformly bounded. This immediately implies that operator f)5.22p is contracting on 
the space of exponentially decreasing functions (<fi(t), p(t))t>o endowed with the norm 
||0, p|| = sup e a *(||v?(t)||^ + k||p(£)||p), for all sufficiently small k, e and 7. This gives us 

the required existence and uniqueness of the fixed point (4>(t), p(£))t>o- Being the fixed 
point of a contracting operator, it depends continuously on every parameter on which 
the operator depends continuously, so it depends continuously on the function V 2 . Note 
that V 2 V 1 as 7 —7- (by (15.201) . this just means that the manifold W y o i£ depends on 
y° continuously). Hence, in the same limit we have (0(t),p(t)) — > (which is the trivial 
fixed point of (J522D when V 2 = V 1 ). □ 

6. Spatially non-walking solutions and their entropy 

In our application to Ginzburg-Landau equation, the p-component of the LDS (14.91) 
describes the temporal evolution of the centers of soliton pairs, namely the deviations of 
the pair centers from the points of a given spatial lattice. This description is valid only if 
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the distances between the soliton pairs are large enough, i.e. the deviations of the soliton 
pairs from the lattice points stay uniformly bounded for arbitrarily large lattice sizes, see 
Section [2j Thus, it is crucial to be able to control the norm of p(t) = {Pk(t)}k=-ooi i- e - 
to keep all pk bounded. 

On the other hand, according to (14.381) . (14. 37p . in the zero order approximation with 
respect to e we have 

(6.1) P*(*)«P*(0) + [ g(y° k (s))ds, 



where y k is a trajectory from the given hyperbolic set A. Thus, an independent diffusive- 
like behavior of the coordinates Pk{t) should be expected |TuVZ07] in the case A is non- 
trivial (chaotic), i.e. the quantities Pk{t) are out of control in this case. 

The main aim of the Section is to show, however, that under some natural assumptions 
on the set A there exists a set of solutions for which at all t € M 

(6-2) l|p(*)ll,<flo 

for some constant Rq » 1. Moreover, this set is large enough, so that it has positive 
space-time entropy. In what follows, in order to simplify notations, we assume that all 
individual ODE's in the uncoupled LDS (14. 8 j) are identical, i.e. /& = / ', g k = g for all k. 



Theorem 6.1. Let the assumptions of Theorem 4-1 hold and let N := dim P. Let us also 
assume that the hyperbolic set A of system (14.21) is transitive and locally-maximal and 
contains N + 1 periodic orbits Z\ : y — Z\{t), Z 2 : y = z 2 {t) ,. . . , Zn + \ : y = z^+iit) with 
periods T\,. . . , Tjv+i respectively. Define the vectors b{ G P , i = 1, . . . , + 1, as follows: 

(6.3) h:=^ [ ' 'g( Zi (t))dt, 

J-i JO 

and require the following properties to be satisfied: 

1. linear combinations of vectors bi generate the whole space P: 

(6.4) P = span{6i,...,6jv-fi}; 

2. there exist strictly positive numbers A{ such that 

(6.5) ASi + A 2 b 2 + ■■■ + A N+1 b N+1 = 0. 

Then, for all sufficiently small e > 0, there exists a uniformly bounded set JC of solutions 
of system (14 ,9p which has strictly positive space-time entropy: 

(6.6) h{K) > 0. 
Proof. We start with the following observation. 

Lemma 6.2. Let (RO) and A6'.5\) hold. Then, for every vector p 6 P, p ^ 0, there exists 
j = J(p) e {1, ■ ■ ■ , N + 1} such that 

(6.7) p ■ b J{p) < 

and, consequently, there exists 5 > such that, for every p ^ 0, 

(6.8) cos(p,6j (p )) < -5. 
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Indeed, suppose there exists p such that (16. 7ft is wrong, i.e. p-bi > for alH = 1, . . . , N+l. 
Multiplying then equality ( 16. 5ft by this p and using that A; > 0, we conclude that p-bi = 
for all i. By (16.41) . this contradicts the assumption p ^ 0. Thus, (16.71) is verified and (16. 8p 
follows immediately from (16 .70 by compactness arguments. 

The idea of the proof of the theorem is as follows. As the set A is transitive and locally- 
maximal, for every two of the periodic orbits Zi and Zj we may choose two different 
heteroclinic orbits Zij m : y = Zij m (t), m — 1,2, that connect them, i.e. 

lim (zi jm (t) - Zi(t + dr )) = 0, lim ( %m (t) - Zj (t + 9±J) = 

t— >— oo J t— >+oo J 

for some constant 9^ m . The orbits Z ijm also belong to A; in fact, the number of different 
heteroclinics is infinite for each pair of periodic orbits in A, but we need only two of them 
for each i and j. The existence of the heteroclinics mean that we may build orbits in A 
which stay for some time near the orbit Zi, then "jump" along any two of the heteroclinics 
Ziji t 2 into a neighborhood of Zj, stay there, then jump again into a neighborhood of 
another periodic orbit, etc.. We will see that for sufficiently small e one can build orbits 
y = {yk(t)}k=-oc °f system (14.91) with a similar behavior for every component yk(t): the 
component stays close to Zi(t) for some time then jumps to Zj(t), etc., moreover the choice 
of the sequence of the periodic orbits the component shadows can be made independently 
for different k. When the component y^ is close to Zi(t) for sufficiently long time, the 
Pfc-component of the associated solution will move in the direction close to bi as time 
grows (see (I6.il) . (l6.3p ). By (16.81) . if the norm of pk becomes large enough we can always 
find a vector bj such that moving in its direction will lead to a decrease in the norm of 
Pf.. Thus, by jumping each time to a properly chosen periodic orbit Zj we may keep the 
norm of all pk bounded. As each jump can be made by at least two different ways (along 
the first or the second heteroclinic) the set of different solutions of system (14.91) we obtain 
in this way will have positive entropy. 

As the first step in implementing this construction we recall the following standard 
result on the "shadowing" in hyperbolic sets. 

Lemma 6.3. There exist 7 > and a > such that for any two orbits y~(t) and y+(t) 
from the hyperbolic set A which satisfy 

(6.9) ||2/_(to)-y+(to)|| <7, 

where 7 > 7 > 0, there exists an orbit y(t) e A and a phase shift 9 such that 
\\y{t)-y-{t)\\ <a f e at fort<t , 

(6.10) 

\\y(t)-y + (t + 9)\\ <C 7 e- Qt fort>t , \9\<C y , 
where C 7 > depends only on 7 and tends to zero as 7 — > 0. 

A proof can be found e.g. in [KaHa95j. The orbit y(t) corresponds simply to the 
intersection of the local unstable manifold of y~(t) with the local stable manifold of y+(t); 
this intersection belongs to A because this set is locally-maximal. 

Combining Lemma 16.31 with Theorem 15.31 we obtain an analogous result for the lattice 
dynamical system (14. 9p . 
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Lemma 6.4. There exist a > and 7 > such that for all e > small enough, for any 
two orbits y± £ A°° such that 

(6.11) ||y _(*o)-y° + (to)l| Y <7, 

where 7 > 7 > and t £ JR., and /or any solution (y_(t), p_(t)) 0/ (I4.9P belonging to the 
invariant manifold W y o e; taere ernst an or&it y° £ A°°, a solution (y(t), p(t)) £ W y o j£ 0/ 
the lattice system (14. 9p . ana 1 tae set 0/ constant phase shifts 9k, k £ Z, snca t/jat 

sup fcez - n fe °+(t + fc )||y < C,e- a ^\ for t > t , 

(6.12) 

sup feez - ^-(*)||y < C 7 e°( t -*») /or t < t , su Pfc6Z ||0 fc || < C7 7 

and 

(6.13) ||y(t) - y_(t)ll* + ||p(t) - p_(t)||, < C ( e a ^\ t < t , 
where C 7 — > +0 as 7 — )■ 0. 

Proof. Indeed, in order to find the required solution (y(t),p(t)), we first construct a 
trajectory y° £ A°°, each component y^if) of which is defined by yjj!_(t) and y° + (t) by 
virtue of Lemma 16.31 such that ( 16 . 1 2[) is satisfied (since the unperturbed system (14. 8 j) is a 
Cartesian product of systems (14. 2p . we only need to apply Lemma [6731 to every component 
in this product). Applying after that Remark 15 A\ we find (in a unique way) the solution 
(y(t),p(t)) of the perturbed system (I4.9p . satisfying (I6.13P for t <t . □ 

We are now ready to complete the proof of the theorem. We will choose sufficiently 
large constants T and R and sufficiently small constants v and fi and construct a sequence 
of sets Ki of solutions of (14. 9p and a sequence of sets K,® of orbits from A°° such that: 

1) for each of the solutions from /Q there exists an otbit y° £ /C° such that the solution 
belongs to the invariant manifold W y o e ; 

2) for every trajectory y° = {yk(t)}k=-oc e f° r ever y k £ Z there are periodic orbits 
Z ik+ : y = z ik+ {t) and Z ik ^ : y = z ik _{t) (from the set of periodic orbits Z 1: ...,Z N+1 
under consideration) such that 

(6.14) \\y° k (lT) - z ik+ (r kl )\\ < v, \\y° k (-lT) - z tk _(r k ^)\\ < v 

for some (irrelevant) constants r £ [0, 7~], where T = max i=lj ..^at+i T (the periods of Z^)\ 

3) at t = ±IT 

(6-15) ||p(t )||p<i2; 

4) for every solution (y(t),p(t)) £ /Q +i there exists a solution (y(t),p(t)) £ /Q such that 

(6-16) ||y(t) - y(t)|| Y + ||p(f) - p(t)|| P < ^ lT ~^ 

for all \t I < IT (the constant a > depends on the set A only). 

By condition 4, the sequence of the sets /Q converges, as I — > +00, to a certain set 
/C of solutions of the LDS (14.91) (convergence is uniform on any bounded time interval). 
Moreover, conditions 2, 3 and 4 imply that solutions in the set K, are uniformly bounded, 
in particular ||p(t)|| is uniformly bounded for all of the solutions. Thus, to prove the 
theorem, we need to actually construct the sequence /Q and to do it in such a way that 
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the sets /Q would contain "sufficiently many" solutions - this would warrant the positivity 
of the space-time entropy of the limit set /C. 

As K,q we choose the set that consists of one orbit y°(t) = {yl = z i(t)}kL-oo> t ne set 
/Co will consist of one solution in the invariant manifold W y o e which satisfies p(0) = 0. 

Now assume we have built the sets /C°, /Q for some /, and let us construct the sets /C° +1 , 
fCi+\. Take any pair {y° G /C z °, (y, p) G K\ H W y o e } . Let i k ± {k G Z) be the sequences 
of indices defined by ( 16.141) and jk± '■= J(p*;(±ZT)), where the integer- valued function 
J(p) is defined by (16. 8ft . Choose any two sequences rrik± {rrik± = 1 or 2). Choose an 
orbit y° G A°° as follows: y° k+ (t) = Zi k+ j k+ m k+ {t - IT + T kl ), where y = z ijm (t) is one of 
the two (chosen above) heteroclinic orbits ^1,2 which connect the periodic orbits Zi and 
Zj. We assume here that the time parametrization on the heteroclinic orbits is chosen 
such that \\zij m {t) — Zi{t)\\ — v at t = T, and ||% m (£) — Zi(t)\\ < v at all t < T. Hence, 
\\zijm( T ki)— Zi(T~ki)\\ < v (recall that the numbers tu are bounded by T) , so ||y°— y+|| Y < 2u 
by ( 16. 14(1 . Therefore, if v is small enough, we may apply Lemma 16.31 (with y° taken as 
the orbit y°_ of the lemma) and obtain a solution (y(t),p(t)) such that 

(6-17) ||y(t) - y(t)\U + || P (t) - p(t)\\ r < 

at t < to = IT; moreover this solution belongs to the invariant manifold W y o £ associated 
with the orbit y° G A°° such that, as t — > +00, the components tend exponentially 
to the heteroclinic orbits Z ik j k+ m k _ hence to the periodic orbits Zj k+ . 

Absolutely analogously (by applying the version of Lemma 16.31 obtained by inversion 
of time) we obtain the existence of a solution (y(£), p(t)) such that 

(6-18) ||y(f) - m\\r + IIpW - P(*)ll, < 

at t > to = —IT; moreover this solution belongs to the manifold W y o associated with 
the orbit y° G A°° such that at each k the component y1{t) tends exponentially to the 
heteroclinic orbit Zj i m as t — > —00 (and it still tends to Zi ~ m as t — > +00). 

By (16.181) . (I6.17p . condition (16. 16j) is fulfilled by the newly constructed solution (y, p). 
Since each component tends to the corresponding periodic orbit Zj k+ as t — > +00 
and to Zj k _ as t — > —00, and the convergence is, by construction, uniform for all k, I and 
for all possible initial solutions (y(t),p(t)) G /Q, it follows that condition ( 16.14(1 will be 
satisfied for the orbit y° at t = ±(Z + 1)T, provided T was chosen large enough. 

It follows, furthermore, that if T is sufficiently large and e is sufficiently small, then the 
change in pk along the orbit (y(t), p(t)) for the time from t = IT to t = (I + 1)T equals to 
TUj where b'j k+ is uniformly close to the vector bj k+ defined by (16.31) . As jk+ = J(pk(lT)) 
and pk(lT) is close to PkilT) (see (16. 16(1 ). it follows that 

cos(p k (lT),b' lk+ )<-5/2, keZ 

(see (16.81) ). Therefore, 

(6.19) P 2 k ((l + 1)T) = (p k (lT)+Tb' jk+ r <pl(/r) + r||4 + ||(T||4 + || -5\\p k \\). 
By (16.161) and (16. 1 5() ||p fe (/T)|| < R + //, and we see now from (16. 19() that 

Pl((l + 1)T) < R\ 
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provided T is taken sufficiently large with respect to jj and R is sufficiently large with 
respect to T (note that \\bj\\ is bounded away from zero by virtue of (I6.4p . (16. 5p ). Analo- 
gously, one checks that 

PU-(i + m<R 2 - 

As we see, condition (16. 15j) is satisfied by the solution (y, p) at to — + 1)T. 

Thus, we have shown that given any solution from the set JCi and any pair of sequences 
(these sequences define which of the two heteroclinic connections is used to jump 
from the periodic orbit Zi k± to Zj k± ) we obtain a solution which satisfies above conditions 
1-3 with I = l + 1, i.e. the newly built solution can be included into the set /Q +i; we 
have also checked condition 4 that ensures the convergence of the sequence of sets /Q as 
/ — > +oo. As we may choose the sequences in an arbitrary way at each step of the 
procedure, the number of solutions in the set /Q which stay at a bounded away from zero 
distance from each other at \t\ < IT and \k\ < n equals to 4'( 2n+1 ). This immediately 
shows that the space-time entropy of the limit set /C is strictly positive. □ 

Note that the assumption that the set A is locally-maximal and transitive can be 
formulated in a more constructive way. Indeed, assume that we have a set of hyperbolic 
periodic orbits Zi, . . . , Zjv+i, which satisfy conditions 1 and 2 of the theorem. Build an 
oriented graph with N + 1 vertexes: the edge connects the vertex i with vertex j if we 
know there exists a heteroclinic orbit Zy which corresponds to a transverse intersection 
of the unstable manifold of Zi with the stable manifold of Zj. If this graph is transitive, 
then the set A of all orbits which stay for all times in a sufficiently small neighborhood of 
the union of the hyperbolic periodic orbits Zi and the transverse heteroclinic orbits is 
uniformly-hyperbolic, transitive and locally- maximal [AfrS74], so Theorem 16. II holds. 

Note also that assumption (I6.5P is really important for the proof of the theorem. Indeed, 
consider the case dimP = 1 for example. Here the integrals b\ and 62 are real numbers, 
and if condition (16. 5p is violated, they both have the same sign, positive, say. In this case, 
when the component y^ stays close to either of the periodic orbits Z 12 , the component 
Pk will increase with time, so we cannot keep Pk{t) bounded by mere switching between 
Z\ and Z 2 . However, assumption (16 .5p can be relaxed if we allow for a uniform drift, 
common for all Pk{t). Namely, the following statement holds true. 

Corollary 6.5. Let all of the assumptions of Theorem \6.1\ be fulfilled except of ( 16. 5p . 
Assume the convex hull lb of vectors bi, ■ ■ ■ b N+ i have a non-empty interior: 

(6.20) I h := int{conv{& 1 , • • • , b N+1 }} ^ 0. 

Let p C h- Then, for every sufficiently small e, there exists a set K,^ of solutions 
(y(t), pit)) of system (14. 9p such that )Cp has positive space-time entropy and each solution 
from /Cp- satisfies 

(6.21) ||y(t)|| ¥ + ||p(t)-pt|| p <i? <oo, teR, 

where the constant Ro depends on p, but is independent oft and the choice of the solution. 

Indeed, for every p E lb conditions (I6.4p .f l63]) hold for the vectors b\ — p, 62 — P, • • • , 
&at + i — p. Then, applying Theorem 16. II to the system obtained from (14. 9 p by substracting 
p from the function g, we immediately obtain the corollary. 
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Remark 6.6. In the one-dimensional case (dimP = 1), we only need two hyperbolic 
periodic orbits, Z + and Z_, connected by transverse heteroclinics. Conditions (16. 4p . (I6.5P 
read now 

(6.22) [ g{z-(t))dt- I + g(z+(t))dt<0; 

Jo Jo 

conditions (16 .4p . (I6.20j) read as 

(6.23) i- £ g{z_{t))dt^±-j\{z + {t))dt. 

In order to establish the existence of the heteroclinic cycle with two hyperbolic periodic 
orbits one may use Shilnikov criterion. Namely, it is enough to show the existence of a 
saddle- focus equilibrium state y = z with a homoclinic loop y = Zh(t), z h {t) — > z as 
t —7- ±oo, and to check that the so-called Shilnikov conditions of chaos are satisfied (we 
will not discuss a higher- dimensional case as in the application we consider in this paper 
we have y G -R 3 ; in the three-dimensional case the Shilnikov condition is that the nearest 
to the imaginary axis characteristic exponent is not real; the equilibrium state must be 
hyperbolic, i.e. it has characteristic exponents on both sides of the imaginary axis and 
no characteristic exponents on the axis). Then there exists a sequence Z n of hyperbolic 
periodic orbits which converge to the homoclinic loop as n — > +oo, any two of them are 
connected by transverse heteroclinics [Shi65llS"hi70| . The periods T m of Z m tend to infinity. 
One can always choose time parametrization such that sup tc ^ r m T m j \z m (t) — Zh(t) — z \ — > 

Oasm-y +oo. It follows that one can always choose among the orbits Z m a pair satisfying 
conition (I6.23j) . provided 

(6.24) / (g(z h (t))-g(z ))dt^0 



(the integral converges since Zh(t) tends to zq exponentially - because of the hyperbolicity 
of zo). Note that the homoclinic loop to a saddle-focus may split as we perturb the system, 
however the two hyperbolic periodic orbits that we find near the loop do not disappear, 
nor the transverse heteroclinics that connect them do, so by checking condition (I6.24p for 
one parameter value we establish the existence of spatio-temporal chaos for an open set 
of parameter values; see Lemma 12.31 in Section [2] for an example. 

Remark 6.7. In the case of LDS with n spatial dimensions, i.e. those parameterized by 
multiindices k G Z n instead of k G Z, the result of Theorem 16.11 obviously, remains true 
under the properly modified definition of the space-time topological entropy. In fact, this 
case is just formally reduced to k G Z by an appropriate reparameterization of the grid 
Z n by the points from Z. 

Appendix. Soliton interaction equations 

We do not reproduce a proof here, however we provide a very informal derivation of the key equations 112,2411 . for the sake 
of completeness. Write equation 1 12.111 in a symbolic form dtu = H(u) +eQ(u), where the operator H is such that H(U) = 0, 
and the small parameters e govern the perturbation, i.e. e = (fi, 5 — So) in our case, and eQ(u) := /i — i(8 — 8o)u. By making 
the multi-pulse ansatz u m ^2 Ug. (t),<f>j (t) ( se e 112. 1911 ). we obtain the following equation: Y2{—e"^^i ipi{x — + e"^> 1^2(3; — 
£i )&j } ~ H(u) + eQ{u) (see 112.411 ). which gives us, after taking the inner product (defined by II2.8I I) with e~ %< ^i ipi t 2 — 
that ij ss - Re /+" e'^itp^x - ij){H + sGjdx, and <j>j « Re e'^i ip 2 (x - £j){H + eG}dx. As the functions ^1,2 
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decay exponentially, and the distances |£j — fi— 1| assumed all to be large, the main contribution to these integrals is given 
by the segment x S + As the pulse U(x) decays exponentially fast, the main contribution to the 

approximate solution u = ^2U^ ky( p k on this segment is given by the term Ug.^. := e"^^U(x and the next order cor- 

rections come from the two terms ±1 ,0 ±1 ■ By recalling that H(U^ ■ =0, we obtain (after shifting x i— ¥ x + £j in the 
integral) & « - Re /+« % (a^i {U (j+1 , 0J + 1 + -J,^ _ + £ Re e"^ ^ (x)g( e ^ C/(aO)dx, 

where £(7 := H' (U) is given by ||2,6[ I. and 0j(x) = 1 at 1 £ [^(Sj— 1 — £j)> § — and vanishes at the other values 

of x. We may rewrite the equation for ^ as £ 3 ~ - Rc/+^ £^{^(3)^ (asJXl/^j-^,^^-^ + U$._ 1 _^ )dx + 

eRe J"+~ e- i<t 'iip 1 (x)g(e i 't'iU(x))dx, where £t is conjugate to C u (see J27TJ). As £t^i = 0, it is easy to see from J2~7t that 
£^ {9j {x)ipi (a;)} = (1 + i/3o)[Sj / (^)'0l (z) + 2£K (a;)i/>j (a)]. By noticing that 6'. (x) is the difference of two delta-functions, the 
one at x = k~ := ^(fj— 1 — £j) minus the other at x = := i (fj+i — one immediately evaluates the first integral in the 
equation for £; to find that the second integral is zero we notice that in our case ee~ '^J Q(e"^i 'U(x)) = /ze — — (S — 5o)tf2 
and use the normalization condition 112,81 1. As the result we get 

Rc(l+il3o)^e t ^3-i-^3^^ 1 (x)U'{x-2K-)-^[(x)U(x-2Kj)] x ^ K - - e < (*i+ 1- *i) [ipi(x)U'(x -2«t)_ ^(a;)Z7(x -2*;+)]^ 
Similarly, 

<^-« - Re{l+ip )^ i ^j-i-^o\i> 2 {x)U\x-2K-)- i>' 2 (x)U{x-2K- )} x=K - - e^i+i-^i) [^2(3)1/' (a; - 2k+)- ^(a;)C/(a; - 2«t)]_ 
+/iRece i( ^- < M - (5 - <S ) 

(see H2.8I I, (??))■ As = |(£j±i — £j) are large in absolute value, we may use asymptotics l|2.10ll . which gives us expressions 
for £j and 0j compatible with H2.22H indeed. Similar derivations can be found e.g. in ITuVZOf] ; a complete proof is in 
IMieZe09l . 
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